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SECTION  I 


INTRODUCTION 


The  Aeroballistic  Research  Facility  Data  Analysis  System  (ARFDAS)  is  a  set  of 
Fortran  programs  used  for  the  determination  of  aerodynamic  coefficients  from 
spark  range  experimental  data.  These  programs  are  intended  for  use  by 
engineers  at  the  United  States  Air  Force  Aeroballistic  Research  Facility 
(AP.F),  an  enclosed  spark  range  used  to  examine  the  ballistic  trajectories  of 
experimental  munitions  in  free  flight  (Reference  1).  Two  of  the  programs. 
Projectile  Design  Analysis  System  (PRODAS)  and  Modified  Trajectory  (MODTRAJ), 
are  designed  to  estimate  the  aerodynamics  of  a  projectile  and  create  simulated 
test  data.  These  programs  are  used  to  "design  the  experiment.”  The  remaining 
programs  use  the  test  data  (actual  or  simulated)  to  determine  the  aerodynamic 
coefficients  from  a  best  fit  of  the  theoretical  equations  to  the  observed 
motion.  The  following  menu  is  a  list  of  program  options  currently  available 
in  ARFDAS. 

Enter  the  number  of  the  desired  action 
-  Linear  Theory  - 

1  -  Input  /  Analysis 

2  -  Output 


6  DOF 


3  -  Setup 
A  -  Analysis 

5  -  Output 

- Auxiliary  — — — - 

6  -  Plot  Experimental  Data  Points 

7  -  Delete  stored  files 

8  -  Run  Time  -  Distance  Bias 

9  -  Run  MODTRAJ 
10  -  Run  PRODAS 

11-  Breakup  Tunnel  XIZ  files 

99  -  Exit 

In  a  typical  analysis,  PRODAS  would  be  run  first.  In  PRODAS,  a  finite  element 
model  of  the  projectile  is  created  from  which  aerodynamic  coefficient 
estimates  of  the  model  are  found.  These  estimates  will  be  used  as  the 
starting  point  for  a  least  squares  fit  of  the  test  data.  An  exterior 
ballistics  analysis  may  be  performed  to  determine  the  projectile  motion.  This 
data  would  be  useful  in  locating  cameras  for  a  test  firing.  PRODAS  may  be 
used  to  create  input  data  for  the  program  MODTRAJ.  In  MODTRAJ  a  Six  Degree  of 
Freedom  (6  DOF)  motion  analysis  is  performed.  The  motion  data  is  correlated 
to  each  of  the  ARF's  camera  locations.  Typical  noise  is  then  superimposed  on 
the  motion  data.  The  resulting  file  may  be  run  through  the  data  reduction 
codes  (linear  and  6  DOF)  to  test  for  coefficient  sensitivity. 
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After  actual  teat  data  has  been  read  from  the  film  the  Time-Distance  Bias 
program  is  run.  This  program  adjusts  the  projectile  motion  to  account  for 
measurement  bias  as  found  from  a  previously  conducted  range  survey 
(calibration). 

The  biased  experimental  data  may  now  be  plotted.  This  visual  examination 
makes  it  easier  to  show  problems  with  the  data  such  as  points  obviously  out  of 
position  or  long  gaps  between  stations. 

The  first  step  in  data  reduction  is  the  linear  theory  input/analysis.  The 
linear  theory  analysis  program  creates  input  for  the  6  DOF  program.  Plotted 
or  tabulated  output  may  be  created  from  either  the  linear  theory  or  6  DOF 
programs.  This  report  will  further  explain  the  theory  of  the  analysis 
techniques  and  Instructions  for  using  the  ARFDAS  programs. 

ARFDAS  is  implemented  on  one  of  the  Eglin  AFB  VAX  computers.  The  program 
utilizes  the  Tektronix  graphics  library  Plot-10. 

Throughout  all  the  programs  that  comprise  ARFDAS,  common  techniques  for 
modifying  data  are  used.  The  principal  method  is  to  display  a  menu  of 
parameters  for  modification  by  the  user.  This  is  done  by  specifying  the  line 
number  in  the  menu  and  the  new  value: 


S87040613  40m  HEDP  Tubular 


1  -  Temperature 

(deg  C): 

0.00 

2  -  Pressure 

(mbars): 

0.00 

3  -  Relative  humidity  (40Z-.40): 

0.00 

4  -  Projectile  mass 

(grams): 

0.00 

5  -  Projectile  diaw 

ster  (inches): 

0.00 

Key  in  item  number 

new  value 

Key  in  "LIST"  to  relist  "DONE"  to  continue 
;  1,19.6 

In  the  above  example,  the  temperature  will  be  changed  from  0  to  19.6  degrees 
Celcius.  The  "LIST"  option  will  redisplay  the  input  menu  showing  any  changes 
that  have  been  made.  The  "DONE"  option  signals  that  modification  is  complete 
on  this  menu.  Only  the  first  character  of  the  command  is  examined,  therefore, 
only  an  L  or  D  must  be  entered.  The  program  will  recognize  either  upper  of 
lower  case  characters. 
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SECTION  II 


SIMULATED  RANGE  DATA  (MODTRAJ) 


MODTRAJ  is  a  modified  6  DOF  trajectory  analysis  program.  Instead  of 
outputting  analysis  results  throughout  the  trajectory  at  specified  time  or 
distance  increments,  the  data  points  are  output  when  the  horizontal  travel 
matches  the  location  of  a  camera  station  in  the  ARF.  Random  noise  is  then 
superimposed  on  the  six  output  parameters  (X,  Y,  Z,  pitch,  yaw,  roll)  before 
the  analysis  results  are  written  to  the  data  file  for  reading  by  Linear  Theory 
Input.  The  probable  errors  of  the  noise  were  determined  from  the  reading 
accuracy  of  the  ARF  range  and  are  hardwired  into  the  code.  The  only  user 
input  into  MODTRAJ  during  program  execution  is  the  name  of  the  data  file 
written  by  PRODAS.  The  output  data  file  from  MODTRAJ  will  have  the  same  name 
as  the  input  file  and  is  in  the  form  of  a  Linear  Theory  input  file. 

The  input  for  MODTRAJ  is  created  by  the  trajectory  segment  of  PRODAS  (Section 
V).  Selecting  the  MODTRAJ  option  in  PRODAS  will  cause  the  contents  of  all 
trajectory  input  menus  to  be  written  into  a  data  file  for  reading  by  MODTRAJ. 
PRODAS  will  not  perform  a  trajectory  analysis  when  running  the  MODTRAJ  input 
option. 

Once  MODTRAJ  has  generated  a  trajectory,  the  resulting  data  points  may  be 
plotted.  The  following  output  menu  shows  the  choices  of  plotted  output: 

Theoretical  Point  Plots 


▼s  time  vs  travel 

1  -  Down  Range  Travel  [X]  -  11 

2  -  Horizontal  Motion  [Y]  -  12 

3  ~  Vertical  Motion  [Z]  -  13 

4  -  Pitch  [THETA]  -  14 

5  -  Taw  [PSI]  -  15 

6  -  Roll  [PHI]  -  16 


7  -  Pitch/Taw  [THETA-PSI] 

8  -  Pitch/Taw  Step  Plot 

9  -  Done 

:  4^ 

Figure  1  shows  a  sample  plot  of  pitch  data  for  a  20mm  projectile.  Tabulated 
data  may  also  be  created  by  MODTRAJ  upon  exit.  Table  1  shows  a  partial 
tabulated  output  created  by  the  20mm  sample  run. 

A  primary  application  of  MODTRAJ  is  to  generate  data  for  simulated  data 
reduction.  This  reduction  can  determine  parameter  sensitivity.  Before 
performing  a  6  DOF  reduction  on  actual  test  data,  the  reduction  of  MODTRAJ 
data  can  determine  which  coefficients  can  be  reduced  with  acceptable  accuracy. 
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20mm  TP  MODTRAJ  test 
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Figure  1.  MODTRAJ  Data  Output 


TABLE  1.  MODTRAJ  TABULAR  OUTPUT 
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SECTION  III 


TIME -DISTANCE  BIAS  (TDBIAS ) 


The  Time-Distance  Bias  (TDBIAS)  program  serves  two  objectives:  (I)  a 
preliminary  fit  to  the  experimental  time-distance  data  to  screen  for  erroneous 
timing  data,  and  (2)  the  option  of  superimposing  statistical  bias  data  based 
on  the  range  calibration.  The  program  will  flag  the  user  to  potentially 
suspect  data.  Corrections  to  the  data  may  be  entered  from  within  the  program. 
Generally  the  statistical  biasing  data  will  be  superimposed  to  improve  data 
quality.  For  biased  data  the  program  will  prefix  the  shot  number  with  the 
letter  B. 

The  data  resulting  from  the  film  readings  is  stored  in  a  file  containing 
multiple  shots.  The  first  step  in  program  execution  will  be  to  specify  the 
name  of  this  file.  The  shots  from  this  file  will  be  selected  individually  for 
biasing.  A  LIST  option  is  available  to  display  the  name  of  all  shots  within 
the  composite  file. 

The  following  illustrates  the  shot  selection  procedure: 

Enter  the  name  of  the  file  containing  the  raw  data:  TUNDAT 
Key  in  10  digit  projectile  title  (or  LIST):  LIST 
Basic  range  data  projectile  titles  (TUNNELXYZ) 

S87040613  40  MM  HEDP  TUBULAR  REREAD 
S87040209  40  MM  HEDP  TUBULAR  REREAD 
S87080617  20  MM  BOOMED /CASED 
S87080618  20  MM  BOOMED/CASED 
S87032490  40  MM  HEDP  TUBULAR 
S87032797  40  1M  HEDP  TUBULAR 

Key  in  10  digit  projectile  title  (or  LIST):  S87040613 

For  the  shot  selected,  the  following  physical  and  atmospheric  properties  must 
be  entered: 

S87040613  40mm  HEDP  Tubular 


1  -  Temperature 

(deg  C): 

0.00 

2  -  Pressure 

(mbars): 

0.00 

3  -  Relative  humidity 

(40Z-.40): 

0.00 

4  -  Projectile  mass 

(grans): 

0.00 

5  -  Projectile  diameter 

(inches): 

0.00 

Key  in  item  number  new  value 

Key  In  "LIST"  to  relist  "DONE"  to  continue 

:  1,19.6 
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A  curvilinear  regression  fit  is  performed  on  the  data  using  the  1MSL 
mathematics  library.  This  curve  fit  finds  a  3rd  order  time  solution.  This 
equation  along  with  the  input  echo,  experimental  times,  and  curve  fit  times 
will  be  printed  out: 

BS87040716  40  MM  HEDP  TUBULAR 

TIME  -  -0.3530054E-10  *  X**3  +  0. 1876861E-06  *  X**2  + 

0.9585853E-03  *  X  +  0. 5469398E+00 

Temperature  (Deg  C):  19.60 

(Deg  F):  67.28 

(Deg  R):  526.92 

Pressure  (abar):  1017.000 
(lbf/ft**2):  2124.513 

Humidity  :  0. 520 

Mass  (grams):  244.00 

Mass  (slugs):  0.016719 

Diameter  (in):  1.575 

Area  (ft**2):  0.01353 

PV  (lbf/ft-sec):  24.74380 

RH0  (slug/ft**3):  0.00233652 

DBSQ  :  0.00 

TA  (deg):  0.00 

BS87040716  40  HEDP  TUBULAR 

Distance  Input  Time  Fit  Time  Delta  Velocity  CD 

(ft)  (sec)  (sec)  (sec)  (ft/sec) 

11.9377003  0.5584267  0.5584099  -0.1680851E-04  1038.366  0.409506 
17.1982994  0.5634965  0.5634812  -0. 1525879E-04  1036.259  0.407454 
27.1723995  0.5731447  0.5731248  -0.1990795E-04  1032.304  0.403592 

47.2350998  0.5925861  0.5926338  0.476837ZE-04  1024.506  0.395938  *** 
57.1082993  0.6022863  0.6022886  0.2324581E-05  1020.745  0.392227 

Standard  Deviation  :  0.0000159308 

XR  Average  :  0.0000000191 

Note  in  the  above  partial  printout  that  the  point  at  47.235  feet  is  flagged 
with  "***."  Any  point  where  the  error  in  time  exceeds  1.5  standard  deviations 
will  be  flagged.  This  flag  is  only  to  call  attention  to  the  point  for 
possible  action  by  the  engineer;  no  program  action  is  taken. 
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Once  Che  curve  fit  has  been  displayed,  Che  user  has  Che  following  opCions: 
Enter  number  of  the  desired  option 

1  ~  Edit  current  data 

2  -  Plot  fit  results 

3  -  Refit  the  current  data 

4  -  Store  the  edited/biased  data 

5  -  Select  a  new  shot 

6  -  Print  tabulated  file  on  exit  or  new  shot 

7  -  Exit 


: 

Editing  may  be  done  on  Che  time  or  roll  at  each  station.  Entering  a  negative 
time  will  remove  the  station  from  the  analysis.  A  value  of  “99  for  roll  will 
result  in  the  roll  data  being  ignored  by  the  Linear  Theory  and  6  DOF  analyses. 
The  following  is  a  partial  listing  resulting  from  entering  the  edit  option: 

Tunnel  XIZ  data  for  editing 


Enter  negative  time  to  remove  a  station 


Station 

Station 

Distance 

Time 

loll 

Distance 

Time 

Roll 

(feet) 

(sec) 

(deg) 

(feet) 

(sec) 

(deg) 

11.9377 

( 

1) 

0.5584267 

(14)  -99. 

242.0256 

(27) 

0.7894410 

(39)  -99 

17.1983 

( 

2) 

0.5634965 

(15)  -99. 

252.0643 

(28) 

0.7999383 

(40)  -99 

27.1724 

( 

3) 

0.5731447 

(16)  -99. 

271.9348 

(29) 

0.8207810 

(41)  -99 

47.2351 

( 

4) 

0.5925861 

(17)  -99. 

276.9530 

(30) 

0.8260718 

(42)  -99, 

57.1083 

( 

5) 

0.6022863 

(18)  -99. 

291.9944 

(31) 

0.8419708 

(43)  -99 

67.0824 

( 

6) 

0.6120735 

(19)  -99. 

306.9884 

(32) 

0.8578804 

(44)  -99, 

72.1105 

( 

7) 

0.6170182 

(20)  -99. 

322.0588 

(33) 

0.8739316 

(45)  -99, 

76.9907 

( 

8) 

0.6218238 

(21)  -99. 

352.1425 

(34) 

0.9062108 

(46)  -99, 

151.9478 

( 

9) 

0.6967998 

(22)  -99. 

382.1923 

(35) 

0.9387400 

(47)  -99 

181.9187 

(10) 

0.7273291 

(23)  -99. 

397. 1584 

(36) 

0.9550332 

(48)  -99. 

186.9524 

(ID 

0.7324781 

(24)  -99. 

412.1839 

(37) 

0.9714601 

(49)  -99, 

212.0999 

(12) 

0.7583820 

(25)  -99. 

442. 1961 

(38) 

1.0044963 

(50)  -99, 

222.1359 

(13) 

0.7687663 

(26)  -99. 

Key  in  item  number  ","  new  value 

Key  in  "LIST"  to  relist  "DONE"  to  continue 


:  13.-1 


In  the  above  example,  the  data  associated  with  number  13  (222.1359  feet)  was 
eliminated. 
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The  following  plots  may  be  created: 

Enter  the  nunber  for  the  desired  operation 

1  -  line  -  fit  plot 

2  -  Velocity  versus  Distance 

3  -  CD  versus  Distance 

4  -  Done 


:  J_ 

The  time  plot  includes  both  time  versus  distance  and  error  at  each  station. 
The  bar  graph  makes  it  easy  to  spot  obvious  errors  such  as  erroneous  input 
times.  See  Figure  2  for  a  sample  TDBIAS  time  fit  plot.  The  other  plots  are 
self-explanatory. 

When  the  data  is  stored  it  will  be  written  in  an  individual  file  by  shot 
number  for  reading  by  the  Linear  Theory  input  program. 


SECTION  IV 


PLOTTING  EXPERIMENTAL  POINTS 


This  program  allows  the  user  to  read  in  the  experimental  data  points,  as  they 
come  from  the  film  readings  or  TOBIAS  program,  and  create  scatter  plots.  The 
points  are  plotted  without  any  reduction  or  analysis  being  performed.  This 
allows  the  user  to  examine  the  data  prior  to  attempting  to  analyze  it.  One 
reason  the  engineer  may  wish  to  do  this  might  be  to  determine  the  degree  of 
angular  motion  present  during  the  flight.  Large  angles  could  alert  the 
engineer  that  non-linear  aerodynamics  might  be  present,  or  very  small  angles 
could  indicate  that  certain  angle  dependent  coefficients  might  be  difficult  to 
determine.  The  following  menu  will  be  displayed  for  plot  choice: 

Linear  Theory  Experimental  Point  Plots 


1  “  Horizontal  Drift  [Y] 

2  -  Vertical  Drop  [Z] 

3  -  Boll  [PHI] 

4  -  Pitch  [THETA] 

5  -  Yaw  [PSI] 


vs  Down-Range  Travel  [X] 
vs  Down-Range  Travel  [X] 
vs  Down-Range  Travel  [X] 
vs  Down-Range  Travel  [X] 
vs  Down-Range  Travel  [X] 


99  -  Done 


Figure  3  shows  the  experimental  yaw  data  for  a  40mm  projectile. 


* 
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Figure  3.  Experimental  Points  Only  Plot 


SECTION  V 


PROJECTILE  DESIGN  ANALYSIS  SYSTEM  (PRODAS) 


1.  GENERAL 

PRODAS,  as  installed  at  ARF  is  a  Fortran  program  used  to  aid  the  engineer 
in  performing  stability  and  exterior  ballistic  analyses.  The  program  uses  a 
common  data  structure  to  pass  information  among  independent  modules  to 
facilitate  rapid  projectile  design  and  analysis.  The  program  uses  interactive 
menus  in  combination  with  Tektronix  Plot-10  graphics  to  provide  a  user 
friendly  environment.  Input/output  may  be  performed  in  either  english  or 
metric  units.  The  program  has  been  used  for  projectiles  from  5. 56mm  to  8 
inches  in  diameter. 

The  first  step  will  be  to  specify  the  terminal  being  used: 

1  -  Large  Screen 

2  -  Snail  Screen 

3  -  Large  Screen  Color  (Tek  4105) 

4  -  Snail  Screen  Shaded  Color  (Tek  4105) 

Choosing  the  large  or  small  screen  will  determine  the  number  of  lines 
printed  on  a  page  and  the  number  of  columns  in  a  line.  The  small  screen  is 
limited  to  26  lines  of  80  columns.  The  large  screen  has  55  lines  of  132 
columns  to  a  page.  In  the  small  screen  mode,  the  text  is  output  as  Tektronix 
character  size  2.  In  the  large  screen  mode  the  text  size  will  be  smaller. 

The  large  screen  mode  will  take  fewer  pages  than  the  small  screen  mode  but  the 
small  screen  mode  will  be  more  readable. 

PRODAS  also  supports  the  Tektronix  4100  series  color  terminals.  Options 
Three  and  Four  will  use  these  terminals  to  make  color  line  plots,  and  use  the 
shaded  areas  capability  to  make  color  shaded  drawings  of  the  projectile.  For 
the  color  shaded  drawing  to  appear  correct,  the  model  must  be  constructed  in  a 
specific  manner.  The  rules  for  model  construction  are  explained  in  the  Edit 
New  (EN)  module. 

After  the  terminal  type  has  been  specified,  the  main  menu  for  the  program 
will  be  displayed.  The  appearance  of  this  menu  is  different  for  the  large  and 
small  screen  modes. 
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Select  code  of  desired  analysis  from  the  following  nenu 

EH  Enter  new  data 
EF  Read  existing  data  file 
EE  Edit  existing  data 

M  Physical  Properties 
S  Stability  Analysis 

T  2/6  DOF  Trajectory 

C  Catalog  Data 
DF  Delete  Existing  File 
B  Exit  PRODAS 

Enter  code  for  desired  operation: 

Each  module  of  PRODAS  is  accessed  by  entering  a  one  or  two  character 
mnemonic*  As  each  module  is  completed,  the  program  will  return  to  this  menu. 

2.  PRODAS  PROGRAM  SEGMENTS 

The  following  paragraphs  detail  each  of  the  program  segments. 

When  running  an  analysis,  data  is  passed  among  the  modules.  This  requires 
that  the  analysis  be  run  in  a  specific  order.  Once  data  has  been  cataloged  it 
will  not  be  necessary  to  run  that  module  again  unless  the  projectile  has  been 
changed. 

The  Physical  Properties  (M)  module  must  be  run  before  any  other  analysis 
module. 


a.  EN  -  Enter  New 

The  first  step  in  running  an  analysis  on  PRODAS  will  be  the 
construction  of  the  projectile  model.  The  projectile  model  is  a  finite 
element  representation  of  the  projectile  to  be  analyzed.  The  model  is 
composed  of  positive  and  negative  density  frustums,  each  called  an  element.  A 
model  may  contain  up  to  ISO  elements.  Most  projectiles  can  be  adequately 
modeled  with  50-75  elements.  The  elements  that  make  up  the  model  are  divided 
into  nine  categories,  each  defined  by  the  physical  function  of  the  element. 
These  categories  are  called  component  codes.  The  component  codes  are: 

-2  Boom 

-1  Sabot 

0  Projectile  Body 

1  Penetrator 

2  Filler 

3  Rotating  Band 

4  Ogive 

5  Boattail 

6  Fin 


When  accepting  a  new  model,  the  program  will  prompt  the  user  for  a 
file  title,  description,  and  input/output  units.  The  program  will  then  accept 
element  data  until  a  <cr>  is  entered  (Table  2).  Each  element  requires  eight 
parameters: 

left  diameter 
right  diameter 
length 

reference  length 
density 
component  code 
radius 
color  code 

All  elements,  except  fins,  are  assumed  to  be  radially  symmetric.  The 
model  must  be  constructed  such  that  the  projectile  points  to  the  right  and  all 
geometry  lies  within  positive  space.  The  reference  length  is  the  distance 
from  the  origin  to  the  left  edge  of  an  element.  It  is  not  necessary  for  the 
base  of  the  projectile  to  be  at  the  zero  reference  length.  When  entering 
data,  it  is  not  necessary  to  enter  trailing  zeros. 

The  color  code  is  used  to  create  color  shaded  models  on  Tektronix  4100 
series  terminals.  The  codes  signal  the  program  as  to  which  patterns  to  use  in 
drawing  the  element.  Color  codes  may  run  from  1  to  16.  If  a  code  is  not 
defined  by  the  user,  one  will  be  assigned  by  the  program  when  the  model  is 
read  in  Enter  File  Module  (EF)  based  on  the  component  code.  For  a  model  to 
look  correct  when  drawn  in  color,  it  is  necessary  to  enter  elements  in  a 
specific  order.  Elements  will  be  displayed  in  numerical  order.  Negative 
density  elements  are  drawn  with  the  background  color  (index  0).  The  negative 
elements  must  be  displayed  before  elements  are  displayed  that  fill  in  this 
volume.  The  easiest  way  to  remember  the  proper  sequence  is  to  create  the 
models  from  the  outside  in.  Create  the  rotating  band  and  sabot  first.  Create 
the  fin,  then  the  body,  penetrator,  and  ogive.  Any  high  explosive  filler  will 
probably  be  last.  Within  the  Edit  Existing  (EE)  module  is  the  ability  to 
renumber  elements  to  correct  errors. 

If  a  fin  element,  component  code  6,  was  entered,  the  program  will 
request  the  fin  details.  These  are  seen  at  the  bottom  of  Table  2.  Under 
normal  operation  these  will  be  on  a  different  screen. 

After  the  fin  has  been  defined,  the  program  will  proceed  to  the  Edit 
Existing  (EE)  module  to  accept  corrections  to  the  model  (Figure  4). 
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TABLE  2.  ELEMENT  DATA 


Enter  9  character  file  title 
:  MAFMODEL 

Enter  descriptive  text 
:  EN  DEMONSTRATION 

Enter  type  of  units  for  file 
1  -  English  2  -  Metric 

:  1 

Enter  element  parameters  seperated  by  commas 

left  dia,  right  dia,  length,  ref.  length,  density,  code,  radius,  color  code 

-2  Boom  ~1  Sabot  0  Body  1  Penetrator  2  Filler  3  Band 

4  Ogive  5  Boattail  6  Fin  7  Tracer 

Letters  may  be  used  for  densities: 

S  Steel,  A  Aluminum,  T  Tungsten,  P  Plastic,  D  DU,  H  HE, 

1:  1,1, 2,0, .1,0, 0,0 
2:  1,0, 2, 2.. 1,4, 0,0 
3:  2, 2, 1,0,. 1,6, 0,0 
4  : 


Enter  fin  group  information 
Fin  Types 

1  -  Rectangular  2  -  Delta 

3  -  Clipped  Delta  4  -  Swept  Rectangular 

Number  of  Fin  Oroups<<4)  :  1 

Enter  element  numbers  in  fin  group  1  separated  by  commas 
:  3 


Fin 

Type 

: 

1 

Pro; 

jectile  Diana te r  at  Fin 

(  xn  )  : 

1.0 

Fin 

Trailing  Edge 

Thickness 

(  in  )  : 

.05 

Fin 

Leading  Edge 

Thickness 

(  in )  : 

.05 

Number  of  Fins  in 

Croup 

: 

4 
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Figure  4.  New  Model-Finner 
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b.  EF  -  Enter  File 


When  the  Enter  Existing  File  module  is  selected,  the  program  will 
prompt  the  user  for  a  projectile  model  file  name.  If  the  name  is  not  known,  a 
list  option  is  available.  After  typing  EF  the  following  instruction  will  be 
displayed: 


Key  in  9  character  projectile  title  (or  "LIST"):  LIST 

Lists  may  be  formed'  in  one  of  two  methods: 

Enter  the  type  of  list  desired 

1  -  Brief 

2  -  Full 

In  the  brief  list  only  the  file  names  are  listed.  In  the  full  list 
the  file  description  is  printed  in  addition  to  the  file  names.  Configurations 
which  have  been  identified  as  inventory  rounds  will  be  accompanied  by  an 
asterisk  (*). 

After  a  file  is  selected  the  units  to  be  used  for  the  input/output  are 
selected.  The  file  contains  a  flag  to  show  the  units  used  at  the  time  the 
file  was  cataloged.  These  units  will  be  the  default  condition: 

Enter  number  for  the  type  of  units  to  use  for  I/O 

1  ~  English 

2  —  Metric 

<cr>  for  default  (english) 

A  drawing  of  the  projectile  (showing  some  key  dimensions)  will  be 
created.  The  program  will  then  accept  a  main  menu  option.  A  carriage  return 
<cr>  will  return  the  program  to  a  display  of  the  main  menu. 

c.  EE  -  Edit  Existing 

The  Edit  Existing  module  allows  the  user  to  make  changes  to  the 
geometric  representation  of  the  projectile. 

Upon  entering  this  module,  a  drawing  of  the  projectile  identifying 
each  element  and  the  editing  menu  will  be  displayed,  (Figure  5).  A  dashed 
line  is  drawn  down  the  axis  of  the  projectile.  Element  numbers  that  appear 
above  this  line  have  positive  densities.  The  element  numbers  that  appear 
below  the  line  have  negative  densities. 

Edit  Option  One  allows  the  user  to  change  the  title  and  description 
for  the  model.  If  a  projectile  has  been  declared  to  be  an  inventory  round,  it 
is  necessary  to  change  the  title  in  order  to  store  any  changes  to  the  model. 

Edit  Option  Two  allows  for  the  addition,  deletion,  or  modification  of 
element  data.  If  this  option  is  selected,  the  user  will  be  prompted  to  enter 
an  element  number. 


Edit  Menu 


1  rile  Title 

2  Elements 

3  Fin  Xnforaation 

4  Seale  Model 

5  Tabulate  Model  Data 

6  Soon 

7  Change  Units 

8  Change  Ref.  Length 

9  Change  Element  * 

10  Change  Ref.  Dia. 

11  Change  Density 

12  Mass  Properties 

99  Return  to  Main  Menu  Fi lenaeet  MAFSPIN 

EE  DEMONSTRATION 


Figure  5.  EE  Option  Menu 
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To  delete  an  element,  enter  the  element  number  as  a  negative.  The 
model  will  then  be  redrawn.  Be  sure  to  examine  the  model  before  making 
additional  changes.  Since  no  gaps  are  left  when  numbering  elements,  some  of 
the  element  numbers  will  now  be  changed. 

To  add  elements  enter  the  word  ADD.  The  program  will  then  display 
element  numbers  tc  be  added  and  accept  new  geometry  until  a  <cr>  is  entered 
with  no  data.  Data  is  accepted  in  the  same  form  as  when  the  model  was  first 
entered;  left  diameter,  right  diameter,  length,  reference  length,  density, 
component  code,  radius,  color  code. 

To  modify  an  element,  enter  that  element  number.  A  menu  will  be 
displayed  showing  the  current  data  (Figure  6). 

Edit  Option  Three  allows  for  modification  of  the  fin  information 
(Figure  7).  Three  entry  menus  are  necessary  to  define  fin  data.  In  each 
case,  the  program  will  display  the  current  data.  A  <cr>  will  leave  the 
current  data  unchanged.  The  first  display  is  the  number  of  fin  groups. 
Remember  that  the  program  is  limited  to  three  fin  groups  but  stability 
analysis  can  only  handle  one  fin  group.  In  the  second  display,  the  elements 
make  up  the  fin  group.  The  program  is  limited  to  five  elements  in  a  fin 
group.  The  third  display  is  the  fin  geometry.  This  data  includes  fin  type, 
thickness  of  leading  and  trailing  edges,  projectile  diameter  at  fin,  and  the 
number  of  fins. 

Edit  Option  Four  performs  a  scaling  of  all  dimensions  in  the  model. 

The  user  will  be  prompted  for  a  scale  factor.  All  lengths  and  diameters  will 
be  multiplied  by  this  factor. 

Edit  Option  Five  creates  a  table  of  the  mass  properties  of  the  current 
model  and  the  element  geometry  of  that  model  (Table  3).  After  the  table  is 
displayed,  a  full  screen  drawing  of  the  projectile  will  be  displayed  (Figure 
8).  If  a  laser  plotter  is  available,  an  option  to  make  a  print  of  this  screen 
will  be  offered  after  a  <cr>  is  typed. 

Edit  Option  Six  allows  the  user  to  zoom  in  on  some  portion  of  the 
screen  to  make  a  magnified  view  of  the  model.  The  terminal's  cursor  is  used 
to  locate  two  points.  The  first  point  will  be  the  center  of  the  new  screen 
window.  The  second  point  will  be  the  right  edge  of  the  screen.  See  Figures  9 
and  10  for  an  example.  Each  terminal  uses  the  cursor  differently  so  check 
terminal  instructions  if  problems  occur.  With  most  terminals  a  point  is 
detected  by  typing  <space  bar>  or  typing  <space  bar>  <cr>.  To  return  to  a 
full  scale  drawing,  type  <Z>  instead  of  <space  bar>. 

Edit  Option  Seven  will  toggle  the  display  units  between  english  and 

metric. 


Edit  Manu 


1  rila  Titla 

2  Claaants 

)  rin  Xnforaation 

4  scala  Modal 

5  Tabulata  Modal  Data 

6  Sooa 
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12  Mass  Propartias 
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Figure  6.  EE  Element  Modification 
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Edit  H«nu 

1  File  Title 

2  Elements 

3  Fin  Information 

4  Scale  Model 

$  Tabulate  Model  Data 
C  Soon 

7  Change  Units 
t  Change  Ref.  Length 
9  Change  Elenent  • 

10  Change  Ref.  Oia. 

11  Change  Density 

12  Mass  Properties 

99  Return  to  Main  Menu 


Enter  nunber  of  fin  groups 
(current  aodel  »  1): 

Enter  elenent  numbers  in  fin  group  1 
(current  aodel  29  I 


1  Fin  Type  2 

2  Proj  Dia  at  Fin  (in):  0.793 

3  Fin  Thickness  (in):  0.0970 

4  Fin  Lead  Edge  Thk  (in):  0.0400 

5  Muaber  of  Fins  4 

Enter  variable  no.  *,*  value 
Enter  "DONE**  to  continue 
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Figure  7.  EE  Fin  Modification 
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Figure  8.  EE  Projectile  Drawing 
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Figure  9*  EE  Zoom  Option 


Edit  Menu 


1  Pile  Title 

2  El*a*nti 

3  Fin  Information 

4  Seal*  Nodal 

5  Tabulate  Model  Data 

6  Soon 

7  Chang*  Units 
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12  Hass  Properties 
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1 

Enter  list  o t  elements 
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Enter  delta  (in):  0.0 


Figure  10.  EE  Reference  Modification 
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Edit  Option  Eight  allows  for  the  modification  of  the  reference  length 
for  a  group  of  elements.  Three  methods  exist  for  defining  the  elements  to 
change  (Figure  10).  The  first  method  is  to  enter  a  list  of  elements  separated 
by  commas.  Element  numbers  will  be  accepted  until  a  <cr>  is  entered  with  no 
elements.  Method  2  accepts  lower  and  upper  numbers  for  a  continuous  range  of 
elements.  Again ,  multiple  lines  of  input  may  be  used,  each  with  a  continuous 
range.  Method  3  accepts  a  range  of  reference  lengths.  Any  element  whose 
current  reference  length  falls  within  this  range  will  be  modified. 

Once  the  elements  have  been  identified,  enter  the  amount  to  change  the 
reference  length.  A  positive  change  will  move  the  model  to  the  right,  a 
negative  to  the  left.  Remember  that  the  model  must  not  extend  into  negative 
space. 


Edit  Option  Nine  is  used  to  swap  the  numbers  of  two  elements.  This  is 
needed  only  to  create  a  shaded  color  drawing  on  a  Tektronix  4100  series 
terminal.  Elements  are  drawn  in  numerical  order.  When  operating  in  the  color 
shaded  mode,  PRODAS  draws  each  element  as  a  color  quadrilateral.  Elements 
with  a  negative  density  are  drawn  with  the  background  color.  All  other 
elements  are  drawn  with  the  color  selected  by  the  user  or  assigned  by  the 
program.  If  the  elements  are  not  ordered  properly,  this  could  cause  an 
element  to  be  displayed  before  the  negative  element  is  displayed  that  makes 
room  for  it,  i.e.,  the  penetrator  being  draw  before  the  sabot.  This  edit 
option  allows  for  a  quick  reordering  of  elements  to  make  a  cleaner  color 
drawing. 


Edit  Option  Ten  allows  the  user  to  assign  a  value  to  the  reference 
diameter  for  the  projectile.  This  will  override  any  value  computed  from  the 
geometry. 


Edit  Option  Eleven  allows  the  user  to  specify  an  old  density  and  a  new 
density.  All  occurrences  of  the  old  density  will  be  changed  to  the  new 
density.  Note  that  positive  and  negative  densities  are  never  considered  to  be 
the  same. 

Edit  Option  Twelve  computes  and  displays  the  mass  properties  of  the 
projectile  as  it  is  currently  defined. 

Edit  Option  Ninety  Nine  returns  the  program  to  the  main  level, 
d.  Other  Auxiliary  Options 

(1)  C  -  Catalog 

This  program  option  will  take  the  data  representing  the 
projectile  model  and  the  analysis  input  parameters  and  store  them  in  a  data 
file  for  use  in  the  future.  The  data  file  will  use  the  file  name  specified 
for  the  current  projectile  (see  EN  and  EE  options).  If  a  file  already  exists 
by  this  name,  the  old  file  will  be  lost.  If  the  inventory  round  flag  was  set 
for  this  model  file  the  program  will  not  allow  the  user  to  catalog  the  file 
without  first  changing  the  name  (see  EE  option).  The  inventory  round  flag  may 
only  be  set/cleared  by  the  user  using  a  system  editor  on  the  projectile  model 
file.  The  flag  may  not  be  changed  from  within  PRObAS. 


28 


When  a  file  is  cataloged,  the  program  will  examine  the  L1ST.MAF 
data  file.  If  an  entry  does  not  exist  for  this  file  name,  an  entry  will  be 
created. 


(2)  B  -  Bye  (exit  PRODAS) 

This  option  is  used  to  terminate  a  PRODAS  session.  The  program 
maintains  flags  for  the  data  that  is  available  for  storage.  If  this  data  has 
not  been  retained,  the  program  will  question  the  user  to  see  if  this  data  is 
to  be  saved  before  exit. 

If  the  projectile  model  is  not  an  inventoried  round  and  an 
analysis  module  has  been  run  that  could  have  changed  stored  parameters,  the 
user  will  be  prompted  to  see  if  a  catalog  option  is  desired  before  exit. 

(3)  DF  -  Delete  File 

This  option  allows  the  user  to  delete  projectile  model  files  that 
have  become  outdated  or  surplused.  The  program  will  prompt  the  user  to 
specify  the  files  to  be  deleted.  Continue  specifying  names  until  finished 
deleting  files,  then  enter  a  <cr>.  If  the  inventory  flag  has  been  set  for  a 
file,  it  may  not  be  deleted.  The  LIST.MAF  file  will  be  updated  to  remove  the 
entry  for  each  file  deleted. 

e.  M  -  Physical  Properties 

The  physical  properties  segment  calculates  the  weight,  center  of 
gravity,  and  moments  of  inertia  for  each  element  in  the  projectile  model. 
Center  of  gravity  locations  are  measured  from  the  reference  length  for  the 
element  (positive  right).  These  individual  properties  are  then  combined  to 
form  the  properties  for  the  entire  composite  projectile.  Within  the  weight 
summaries,  total  refers  to  the  sum  of  all  elements  (interior  ballistic  weight) 
and  projectile  refers  to  the  flying  portion  of  the  model  (exterior  ballistic 
weight).  Inertias  do  not  include  sabot  elements;  therefore,  it  is  an  exterior 
ballistic  property.  Depending  on  the  projectile  configuration,  the  sabot, 
penetrator,  and  HE1  weights  may  be  computed  and  tabulated  separately.  If  fins 
are  present  on  this  model,  the  fin  group  information  is  displayed.  Tables  4 
and  3  show  typical  results  for  a  spin  and  fin  stabilized  projectile. 

The  only  input  required  for  this  segment  is  the  projectile  model  data. 
The  output  from  this  segment  is  used  by  most  of  the  other  PRODAS  modules. 

Once  the  physical  properties  have  been  calculated,  the  user  has  the 
option  of  overriding  the  computed  properties  for  individual  elements  or  groups 
of  elements.  This  is  done  by  selecting  the  "MOD"  option  after  the  output  has 
been  displayed.  If  the  "MOD"  option  is  selected,  enter  the  element  number(s) 
to  be  changed.  Then  enter  the  weight,  inertias,  reference  length,  and  center 
of  gravity  for  the  group  of  elements.  These  properties  will  be  displayed  as 
the  first  element  given  with  the  properties  for  the  remaining  elements  zeroed. 
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TABLE  4.  SPIN  STABILIZED  MOMENTS  (CONCLUDED) 
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TABLE  5.  FIN  STABILIZED  MOMENTS  (CONCLUDED) 
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Boom  -1  -  Sabot  0  -  Projactila  Bjdy 

Panatrator  2  -  HI  3  -  Band 

Ofliva  5  -  Boattail  6  -  T ins 

Tracer 


It  must  be  noted  that  the  weights  and  inertias  given  for  the  composite 
projectile  are  for  the  "flying"  portion  only.  Sabot  elements  and  rotating 
band  elements  on  a  saboted  projectile  are  removed  before  making  these 
calculations.  The  total  weight  is  carried  separately  for  use  in  the  internal 
ballistics  segment. 

f.  S  -  Stability 

The  stability  segment  is  broken  into  three  phases: 

1  -  Computation  of  aerodynamic  and  stability  parameters 

2  -  Empirical  interior  ballistics 

3  -  Plotted/tabulated  output 

The  form  of  the  input  and  output  menus  will  be  different  for  spin 
and  fin  stabilized  projectiles.  Physical  property  data  must  be  available  from 
the  "M"  segment.  The  input  menu  for  a  spin  stabilized  projectile  is: 

Stability  Inputs 


1  ~  Projectile  Diameter 

(inch): 

0.531 

2  -  Meplat  Diameter 

(inch): 

0.050 

3  - 

(calibers): 

0.094 

4  -  Rotating  Rand  Diameter 

(inch): 

0.531 

5  — 

(calibers): 

1.000 

6  -  Ogive  Radius 

(inch): 

531.000 

7  - 

(calibers): 

1000.000 

8  -  Boom  Length 

(inch): 

0.000 

9  - 

(calibers): 

0.000 

10  -  Projectile  Length 

(inch): 

2.800 

11  - 

(calibers): 

5.273 

12  ~  Ogive  Length 

(inch): 

1.550 

13  - 

(calibers): 

2.919 

14  -  Boattail  Length 

(inch): 

0.000 

15  - 

(calibers): 

0.000 

16  —  CC  from  Rose 

(inch): 

1.890 

17  - 

(calibers): 

3.560 

18  -  Rifling  Twist 

(inch/rev): 

23.504 

19  -  (calibers/rev): 

23.862 

20  -  Gun  Bore  Diameter 

(inch): 

0.985 

21  -  Total  Weight 

(lbm): 

0.230 

22  -  Axial  Moment 

(lbm~in**2): 

0.00758 

23  -  Transverse  Moment 

(lbm-in**2): 

0.05857 

24  -  Air  Temperature 

(Deg  F): 

59.000 

Key  in  item  number  new  value 

Key  in  "LIST"  to  relist  "DONE"  to  continue 

Key  in  "PRINT"  for  a  Line  printer  output 
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The  input  menu  for  a  fin  stabilized  projectile  is: 

Finned  Projectile  Aerodynamic  Coefficient  and 
Stability  Determination 


1  -  Number  of  Fins 

• 

• 

6 

2  -  Fin  Type  (1  Rectangular 

2  Delta)  : 

2 

(3  Clipped 

4  Swept) 

3  -  Fin  Soot  Chord 

(calibers): 

2.867 

4  -  Fin  Tip  Chord 

(calibers): 

0.000 

5  -  Fin  Height 

(calibers): 

1.078 

6  -  Fin  Thickness 

(calibers): 

0.070 

7  —  Fin  Lead  Thickness 

(calibers): 

0.057 

8  -  Boattail  Length 

(calibers): 

2.867 

9  -  Boom  Length 

(calibers): 

0.000 

10  -  Boom  (or  Aft  Boattail)  Dia  (calibers): 

1.000 

11  -  Nose  Radius 

(calibers): 

30.535 

12  -  Me plat  Diameter 

(calibers): 

0.043 

13  -  Hose  Length 

(calibers): 

3.308 

14  -  Hose  to  CG  Distance 

(calibers): 

7.204 

15  -  Projectile  Diameter 

(inch): 

1.395 

16  -  Projectile  Cylindrical  Length  (cal): 

7.505 

17  -  Projectile  Base  to  Fin  Trailing  : 

0.000 

Edge  Distance  (calibers) 

18  ~  Axial  Inertia 

(lbw-in**2): 

1.64276 

19  -  Transverse  Inertia 

(lbm-in**2): 

130.69066 

20  -  Height 

(lbm): 

8.206 

21  -  Projectile  Velocity 

(ft/sec): 

5500.000 

22  -  Magnus  Moment  Coefficient 

• 

• 

5.000 

23  -  Roll  Induced  Side  Moment 

Coefficient  : 

0.000 

24  -  Temperature 

(Deg  F): 

59.000 

25  -  Gun  Bore  Diameter 

(inches): 

4.  128 

26  —  Any  T  Pads  on  Fins  (0 

-  no  1  -  yes): 

0. 

Key  in  item  number  new  value 

Key  in  "LIST”  to  relist  "DONE"  to  continue 

Key  in  "PRINT"  for  a  Line  printer  output 


'  Note  the  SPINNER  data  may  be  entered  in  either  length  units  or 
calibers  but  FINNER  is  limited  to  calibers  only.  A  default  ogive  radius  of 
100  calibers  is  used  if  the  ogive  radius  is  set  to  zero.  SPINNER  requires 
that  an  ogive  exists  in  order  to  run  the  stability  analysis. 

Within  the  "T"  segment  the  aerodynamic  coefficients  may  be  modified  to 
reflect  any  experimental  data  that  becomes  available.  When  this  option  is 
exercised  a  flag  is  set  in  PRODAS.  If  that  flag  is  set  when  the  stability 
segment  is  run,  the  user  will  be  prompted  with  the  message: 

This  round  has  previously  modified  aerodynamic  coefficients. 

Do  you  wish  to  recalculate  the  aerodynamic  coefficients  or 
proceed  directly  to  the  stability  analysis? 

Enter  0  to  recalculate  aeros 

Enter  1  to  proceed  directly  to  stability 
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If  Option  One  is  chosen,  the  values  from  the  data  base  will  be  echoed 
in  the  output.  If  option  zero  is  selected,  or  the  empirical  data  has  not  been 
modified,  a  display  of  tlu  flying  (no  sabot  elements)  portion  of  the  pro¬ 
jectile  model  together  with  the  stability  analysis  model  will  be  displayed. 

The  stability  analysis  model  is  generated  from  the  geometry  as  defined  in  the 
stability  input  menu.  This  drawing  is  displayed  for  the  user  to  examine  for 
model  verification  before  running  the  analysis.  Figure  11  shows  the  veri¬ 
fication  model  for  a  typical  FINNER  analysis.  In  this  figure  both  the 
projectile  as  modeled  and  the  aerodynamic  representation  of  that  projectile 
are  displayed.  If  the  analysis  model  does  not  match  the  desired  configura¬ 
tion,  return  to  the  input  to  correct  the  model.  Once  the  user  is  satisfied 
with  the  model,  the  analysis  will  be  run. 

For  both  fin  and  spin  stabilized  projectiles,  the  tabulated  output 
will  include  an  echo  of  the  input  values,  the  table  of  aerodynamic 
coefficients,  and  the  stability  analysis  parameters.  For  a  spin  stabilized 
projectile,  a  table  of  McDRAG  coefficients  will  also  be  printed.  Table  6 
shows  the  output  for  a  spin  stabilized  projectile.  The  parameters  output  in 
this  table  are: 

SPINNER  aerodynamic  coefficients  are: 

Mach  -  Mach  number 

CX  -  Zero  yaw  drag  coefficient 

2 

CX2  -  Yaw  drag  coefficient  per  sin 

CNa  -  Normal  force  coefficient  per  sin 

Cma  -  Pitching  moment  coefficient  per  sin 

CPN  -  Normal  force  center  of  pressure  (calibers  from  nose) 

Cypa  -  Magnus  force  coefficient  derivative  per  sin 

Cnpa  -  Zero  yaw  magnus  moment  coefficient  derivative  per  sin 

3  — 

Cnpa3-  Cubic  magnus  moment  coefficient  derivative  per  sin 

Cnpa5-  Quintic  magnus  moment  coefficient  derivative  per  sin^ 

CPF1  -  Center  of  pressure  of  Magnus  force  at  1  degree  of  yaw 
(calibers  from  nose) 

CPF5  -  Center  of  pressure  of  Magnus  force  at  5  degrees  of  yaw 
(calibers  from  nose) 

Cnpa-5-  5  degree  secant  slope  of  magnus  moment  coefficient 
derivative  (@  5  deg  yaw)  per  sin  a 

Cmq  -  Pitch  damping  coefficient 

Clp  -  Spin  deceleration  coefficient 
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Figure  11*  Geometry  Validetion  Model 
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TABLE  6.  SPINNER  OUTPUT  (CONCLUDED) 
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McDRAG  aerodynamic  coefficients  are; 

MACH  -  Mach  number 

CDO  -  McDRAG  computed  zero  yaw  drag  coefficient 

CX  -  SPINNER  computed  zero  yaw  drag  coefficient 

DELTA-  Difference  between  McDRAG  and  SPINNER 

CDH  -  Transonic  pressure  head  drag 

CDSF  -  Skin  friction  drag 

CDBND-  Rotating  band  drag 

CDBT  -  Boattail  drag 

CDB  -  Base  drag 

SPINNER  stability  analysis  parameters  are: 

MACH  -  Mach  number 

GYRO  -  Gyro  stability  factor 

SBAR  -  Dynamic  stability  factor  @  1  degree  of  yaw 
RECIP-  Dynamic  reciprocal  factor  @  1  degree  of  yaw 
SBAR5-  Dynamic  stability  factor  @  5  degrees  of  yaw 
RECIP5-  Dynamic  reciprocal  factor  @  5  degrees  of  yaw 
SPIN  -  Spin  rate  (rad/sec) 

W1  -  Nutation  frequency  (rad/sec) 

W2  -  Precession  frequency  (rad/sec) 

LI  -  Nutation  damping  factor  @  1  degree  of  yaw  (1/ft) 

L2  -  Precession  damping  factor  @  1  degree  of  yaw  (1/ft) 

Ll-5  -  Nutation  damping  factor  @  5  degrees  of  yaw  (1/ft) 

L2-5  -  Precession  damping  factor  @  5  degrees  of  yaw  (1/ft) 

DELT  -  Integration  time  step  (sec) 

DISP  -  Dispersion  factor  per  .001  inches  of  barrel 
clearance  (mils) 
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Table  7  shows  the  tabulated  output  from  a  F1NNER  analysis.  The 
parameters  output  are: 

FINNER  aerodynamic  coefficients  are: 

Mach  -  Mach  number 

CD  -  Zero  yaw  drag  coefficient 

CNa  -  Normal  force  coefficient  per  sin  CL 

CPN  -  Normal  force  center  of  pressure  (calibers  from  nose) 

Cma  -  Pitching  moment  coefficient  per  sin  CL 

Cmq  -  Pitch  damping  coefficient 

CXB  -  Body  alone  zero  yaw  axial  force  coefficient 

CXF  -  Fin  alone  zero  yaw  axial  force  coefficient 

CnaB  -  Body  alone  normal  force  coefficient  per  sin  a 

CnaF  -  Fin  alone  normal  force  coefficient  per  sin  CL 

CPNB  -  Body  alone  normal  force  center  of  pressure 
(calibers  from  nose) 

CPNF  -  Fin  alone  normal  force  center  of  pressure 
(calibers  from  nose) 

FINNER  stability  analysis  parameters  are: 

P  -  Spin  rate  (rad/sec) 

GYRO  -  Gyro  stability  factor 

SIGMA-  Inverse  of  the  gyro  stability  factor 

RATE  -20/  (l-<7)  rate  of  nutation  wrt  precession  divided  by 
the  rate  of  precession 

P/Wl  -  Spin  divided  by  the  nutation  frequency 
PD/2V-  Non-dimensional  spin  rate 
W1  -  Nutation  frequency  (rad/sec) 

W2  -  Precession  frequency  (rad/sec) 

LI  -  Nutation  damping  factor  (1/ft) 
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TABLE  7.  FINNER  OUTPUT  (CONCLUDED) 
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FINNER  stability  analysis  parameters  -  (continued) 

L2  -  Precession  damping  factor  (1/ft) 

AF  -  Amplification  factor 

YR  -  Yaw  of  repose  (deg) 

Additional  outputs  on  this  page  include: 

For  rectangular  fines: 

AR-  Fin  aspect  ratio 

LAMBDA-  Body  diameter  divided  by  total  fin  span 
(tip  to  tip) 


For  delta,  clipped  delta,  or  swept  fins: 

Epsilon  -  Angle  from  horizontal  to  a  line  between  forward 
point  of  root  chord  to  aft  point  of  tip  chord 

LAMBDA  -  Body  diameter  divided  by  total  fin  span 
(tip  to  tip) 

Once  the  aerodynamic  coefficients  and  parameters  have  been  printed, 
the  user  will  be  asked  to  select  a  muzzle  velocity  and  aircraft  velocity: 

1  -  Muzzle  Velocity  (ft/sec):  4350,000 

2  -  Aircraft  Velocity (ft /sec):  0.000 

Key  in  item  number  new  value 

Key  in  "CALC"  to  calculate 

Key  in  "LIST"  to  relist  "DOME"  to  continue 

Key  in  "PRINT"  for  a  Line  printer  output 


The  aircraft  velocity  may  be  zero  but  the  muzzle  velocity  must  be 
greater  than  zero.  Since  the  muzzle  velocity  is  used  to  compute  the 
projectile  spin,  the  velocity  components  must  be  dealt  with  separately. 


The  user  may  select  the  option  "CALC"  which  will  cause  PR0DAS  to  branch  into 
the  Empirical  Interior  Ballistic  program.  The  Empirical  Interior  Ballistics 
input  menu  is: 


1  -  Peak  Pressure 

2  -  Barrel  Length 

3  -  Propellant  Weight 
or 

4  -  Load  Density 

5  —  Chamber  Volume 


(psi):  60000.0000 

(in):  64.0000 

(grains):  1569.4000 


(lbm/in**3): 

(in**3): 


0.0380 

5.9000 


Key  in  item  number  new  value 

Key  in  "LIST"  to  relist  "DONE"  to  continue 

Key  in  "PRINT"  for  a  Line  printer  output 
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The  program  will  take  the  values  for  peak  pressure,  barrel  length, 
propellant  weights,  and  chamber  volume  to  estimate  a  muzzle  velocity.  The 
results  of  the  analysis  will  be  displayed  in  tabulated  form: 

Interior  Ballistics  Results 

Gun  Parameters 

Barrel  Length  (iu):  64.00 

Chamber  Volume  (in**3):  5.900 

Bore  Area  (in**2):  0.7849 

Projected  Parameters 

Projected  Weight:  0.297  (lbm)  2084.  (grains)  135.0  (grams) 
Projected  Diameter  (in)  0.985 

Ballistic  Parameters 
Charge/mass  Ratio  0.7528 

Expansion  Ratio  9.5139 

Propellant  Weight  -  0.224(lbm)  1569. (grains)  101.65(grans) 
Loading  Density(lbm/in**3)  0.0380 

Ihizzle  Velocity  (ft /sec)  4375.2 

If  the  user  selects  to  send  the  input  menu  to  the  line  printer  file, 
the  output  will  also  be  sent  to  the  line  printer  file.  The  program  will  then 
return  to  the  point  where  the  muzzle  velocity  input  was  displayed.  Note  that 
the  muzzle  velocity  has  not  been  changed  by  the  Empirical  Interior  Ballistic 
program.  If  the  user  wants  to  use  the  value  computed,  it  must  be  entered  at 
this  point. 

Once  the  velocities  have  been  established  a  summary  of  the  stability 
results,  together  with  a  drawing  of  the  flying  portion  of  the  projectile 
model,  will  be  displayed.  Additional  output  options  are  available  at  this 
point  (Figure  12). 

Option  One  will  cause  a  return  to  the  stability  input  menu. 

Option  Two  will  cause  a  return  to  the  muzzle  and  aircraft  velocity 
input  menu. 

Option  Three  will  display  the  menu  of  plots  available  for  output. 

The  menu  available  depends  on  the  form  of  projectile  being  analyzed.  For  fin 
stabilized: 


Stability  Plot  Segment 

vs  Spin/Nut.  Freq. 

11  W1 

12  W2 

13  LI 

14  L2 

15  AF 

20  Return 


vs  Mach  Number 

1  CZ 

2  CNa 

3  CPN 

4  Caa 

5  Cmq 


45. 
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Figure  12.  Stability  Summary 


For  spin  stabilized: 

Stability  Plot  Segment 

1  CX  8  CPN 

2  CNa  9  Gyro  Stability  Factor 

3  Cma  10  Dynamic  Stability  Factor  (SMR) 

4  Cmq  11  LI 

5  Clp  12  L2 

6  Cnpa 

20  Return 

21  Dynamic  Stability  Criteria 


Figure  13  shows  the  normal  force  (CNa)  curve  found  by  this  analysis. 
The  Dynamic  Stability  Criteria  plot  (21)  shows  the  projected  stability/ 
instability  for  each  of  the  Mach  numbers  analyzed  and  where  they  fall  relative 
to  the  stability  curve  (Figure  14).  In  this  example  we  see  that  this 
projectile  will  be  unstable  at  subsonic  speeds. 

Option  Four  will  display  the  plots  available  for  cross  plotting.  In 
stability,  unlike  the  other  PRODAS  segments,  cross  plotting  will  plot  the 
current  analysis  against  data  stored  in  previously  cataloged  analyses  and  not 
against  analyses  made  during  this  session.  If  the  user  selects  to  enter 
stability  without  first  entering  model  geometry,  the  program  will  branch  to 
the  stability  cross  plot  module.  The  data  available  for  cross  plotting  (both 
spin  and  fin  stabilized)  is: 

1  -  CX 

2  -  CNa 

3  -  Caa 

4  -  Cmq 

5  -  Clp 

6  -  Return 

Figure  15  shows  a  cross  plot  of  drag  data  for  spin  stabilized 
projectiles. 

Option  Five  will  return  the  program  to  the  main  PRODAS  menu. 

Option  Six  will  cause  the  printed  output  that  has  been  displayed  on 
the  screen  to  be  written  into  a  file  for  printing  on  a  line  printer. 

Option  Seven  is  available  only  if  a  laser  printer  is  available  that 
supports  PLOT-10.  This  will  cause  a  plot  of  this  stability  summary  to  be 
made. 

If  the  user  enters  the  stability  option  without  entering  a  projectile, 
the  program  assumes  that  the  user  wishes  to  create  cross  plots  from  the 
cataloged  data. 
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Figure  13.  Stability  Output  Plot 
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g.  T  -  2  /  6  Degree  of  Freedom  Trajectory 

This  module  performs  either  a  2  or  6  DOF  trajectory  analysis.  If  the 
module  is  selected  without  first  defining  a  projectile  model,  the  program  will 
proceed  to  the  cross  plot  segment. 

Three  input  menus  are  needed  to  specify  the  basic  trajectory  physical 
and  analysis  parameters.  The  first  menu  is  used  to  define  the  overall 
projectile  configuration,  which  determines  the  additional  input  menus  needed, 
and  what  additional  output  is  desired: 

Trajectory  Analysis  Options  Selection 


1  -  Type  of  projectile  :  2 

(1  -  Fin  stabilised  2  —  Spin  stabilized) 

2  -  Type  of  trajectory  slnulation  2 

(2  -  Two  DOF  6  -  Six  DOF) 

3  -  Perform  FUMER  analysis  0 

(0  -  No  1  -  Tea) 

4  -  Does  the  trajectory  originate  at  an  aircraft:  0 

(0  -  No  1  -  Yes) 

5  -  Perform  internal  ball  rotor  simulation  :  0 

(0  -  No  1  -  Yes) 

6  -  Enter  number  of  barrels  in  gun  1 

7  -  Examine  or  modify  the  aerodynamics  0 

(0  -  No  1  -  Yea) 

8  -  Generate  cross  plot  file  :  0 

(0  -  No  1  -  Yes) 

9  -  Write  a  QUICK6  or  M0DTRAJ  datafile  :  0 

(O  -  No  file  1  -  QUICK6  2  -  MODTRAJ) 


Key  in  item  number  new  value 

Key  in  "LIST"  to  relist  "DONE"  to  continue 

Key  in  "PRINT"  for  a  Line  printer  output 

Line  one  specifies  the  projectile  stability  configuration,  fin  or 
spin.  The  program  will  examine  the  projectile  model  for  fin  elements  to 
initialize  this  parameter. 

Line  two  specifies  which  equations  of  motion  are  to  be  used  for  the 
analysis,  2  or  6  DOF.  By  default  the  program  will  perform  a  2  DOF  analysis. 
For  the  2  DOF  analysis,  aerodynamic  jump,  velocity  jump,  and  yaw  of  repose 
effects  are  taken  into  account.  Aerodynamic  and  velocity  jump  effects  are  the 
result  of  firing  around  an  axis  that  does  not  match  the  direction  of  air  flow, 
i.e.,  side  firing  from  an  aircraft.  This  causes  an  angle  of  attack  for  the 
projectile  greater  than  it  would  normally  have  for  no  yaw  motion.  Yaw  of 
repose  is  the  resistance  of  a  spin  stabilized  projectile's  nose  to  following 
the  ballistic  trajectory,  i.e.  keeping  a  nose  up  orientation  as  the  projectile 
comes  over  the  apex  of  the  arc  in  the  trajectory. 

If  tracer  burn  or  rocket  effects  are  to  be  included,  line  three 
(fumer)  should  be  flagged. 
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For  a  projectile  launched  from  an  aircraft,  line  four  will  cause  the 
aircraft  parameter  menu  to  be  displayed. 

Line  five  will  cause  the  ball  rotor  fuze  input  menu  to  be  displayed. 
This  enables  the  effect  of  a  moving  component  in  the  fuze  to  be  included  in 
the  equations  of  motion. 

The  number  of  barrels  is  specified  in  line  six.  A  multiple  barreled 
gun  is  assumed  to  be  a  gatling  gun. 

The  Stability  analysis  segment  of  PRODAS  computes  aerodynamic 
coefficients.  Line  seven  allows  the  user  to  display  and  modify  these 
parameters. 

Line  eight  allows  the  user  to  flag  the  program  to  store  trajectory 
data  for  the  creation  of  multiple  analysis  plots.  This  trajectory  data  will 
also  be  stored  for  use  by  other  analysis  modules  within  PRODAS  (i.e..  Target 
Penetration). 

Line  nine  allows  the  user  to  specify  that  the  trajectory  input  data  is 
to  be  stored  for  use  by  either  the  QUICK6  or  MODTRAJ  programs.  In  either  of 
these  cases,  trajectory  analysis  will  not  be  performed.  QUICK6  is  a  batch  6 
DOF  analysis  program.  MODTRAJ  is  a  program  for  simulating  a  ballistic  range. 
The  program  creates  data  points  to  simulate  the  camera  X-Y-Z-theta-psi  data 
with  noise.  (See  Section  II.) 

The  menus  flagged  by  the  above  options  will  be  explained  in  the 
following  pages. 

The  second  and  third  menus  specify  the  projectile  physical  properties 
and  basic  trajectory  analysis  parameters: 

Trajectory  Analysis  Conditions  (Page  1) 


1  -  Projectile  Diameter  (in):  0.S31 

2  -  Axial  Inertia  (lbm-in**2):  0.007581 

3  —  Transverse  Inertia  (lba-in**2):  0.058568 

A  -  Center  of  Gravity  (calibers  from  nose):  3.560 

5  -  Drag  Porn  Factor  1.000 

6  -  Fin  Cant  Angle  (deg):  0.0000 

7  -  Number  of  Fins  0. 

8  -  Projectile  Weight  (lbn):  0.2305 

9  -  Switch  time  from  6  DOF  to  2  DOF  999.000 

10  -  Linear  theory  stability  analysis  0 

(0  -  no  1  -  yes) 

11  -  Aerodynamic  jump  included  0 

(0  -  no  1  -  yes) 

12  -  Limiting  total  angle  of  attack  for  which 

nonlinear  coef's  are  considered  (deg):  10.000 

13  -  Gun  Elevation  Angle  (deg  +*up):  1.000 

(0.0  if  aircraft  originated) 

14  -  Gun  Azimuth  Angle  (deg  +"left):  0.000 
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Trajectory  Analysis  Conditions  (Page  1)  -  (Continued) 


15  -  Missile  angle  of  attack-THETAM  :  0.00 

(deg  +*nose  down) 

16  -  Missile  angle  of  sideslip-PSIM  :  0.00 

(deg  +"nose  left) 

Key  in  iten  nuaber  new  value 

Key  in  "LIST”  to  relist  "DONE"  to  continue 

Key  in  "CALC"  to  calculate 

Key  in  "PRINT"  for  a  Line  printer  output 

Lines  one  through  nine  specify  the  physical  parameters.  These  values 
are  calculated  in  the  Physical  Properties  segment. 

Line  ten  allows  for  calculation  of  the  linear  theory  gyroscopic  and 
dynamic  stability  factors. 

Line  eleven  causes  the  aerodynamic  jump  effects  to  be  included. 

Line  twelve  specifies  the  upper  limit  of  angle  of  attack  for 
non-linear  coefficients.  By  default  this  value  is  10  degrees. 

Lines  thirteen  and  fourteen  define  the  gun  orientation. 

Lines  fifteen  and  sixteen  define  the  angular  orientation  of  the 
projectile  at  muzzle  exit.  These  parameters  do  not  have  an  effect  in  a  2  DOF 
analysis.  The  "CALC"  option  will  run  the  Muzzle  Exit  segment  of  PR0DAS  to 
assist  in  the  selection  of  initial  pitch  and  yaw  values.  After  running  Muzzle 
Exit,  the  program  will  return  to  this  input  menu: 

Trajectory  Analysis  Conditions  (Page  2) 


* 


1  -  Pitch  Angular  Velocity 

• 

• 

0.00 

(rad/sec,  +~nose  down) 

2  -  Taw  Angular  Velocity 

m 

• 

0.00 

(rad/sec,  +™nose  left) 

3  -  Initial  Axial  Spin 

• 

• 

13947.201 

(rad/sec,  +*right  hand 

twist) 

4  -  At  nos.  Sea  Level  Tenp. 

(Deg  F): 

59.00 

5  —  Hind  Range  (fps. 

♦■down  range): 

0.0000 

6  -  Cross  Hind 

(fps,  +*left): 

0.0000 

7  -  Vertical  Hind 

(fps,  +-up): 

0.0000 

8  -  Muzzle  Velocity 

(ft/sec): 

4350.0000 

9  -  Initial  Altitude 

(ft): 

0.0000 

10  -  Initial  Horizontal  Range  (ft): 

64.4461 

11  -  Initial  Cross  Range 

(ft,  +*left): 

0.0000 

12  -  Initial  Tine 

(sec): 

0.0000 

13  -  Integration  Tine  Increi 

sent  (sec): 

0.0200000 

if  0.0  default  ■  nutation  freq./20. 

14  -  Maxi  non  Tine 

(sec): 

70.0000 

15  -  Maxinun  Slant  Range 

(ft): 

10000.00 
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Trajectory  Analysis  Conditions  (Page  2)  -  (Continued) 


Output  Increment  Time  or  Slant  Range 

16  -  Time  (sec):  0.1000 

17  -  Slant  Range  (ft):  0.00 

Key  in  item  number  ,,>M  new  value 

Key  in  "LIST"  to  relist  "DONE"  to  continue 

Key  in  "CALC"  to  calculate 

Key  in  "PRINT"  for  a  Line  printer  output 

Lines  one  and  two  define  the  angular  velocity  of  the  projectile  at  the 
nuzzle  exit.  These  parameters  do  not  have  an  effect  in  a  2  DOF  analysis.  The 
"CALC"  option  will  run  the  Muzzle  Exit  segment  of  PR0DAS  to  assist  in  the 
selection  of  initial  pitch  and  yaw  rate  values.  After  running  Muzzle  Exit, 
the  program  will  return  to  this  input  menu. 

Line  three  defines  the  projectile  spin.  A  value  of  60  rad/sec  is 
assigned  for  fin  stabilized  projectiles.  The  program  cannot  handle  a 
projectile  with  a  spin  of  zero.  The  spin  is  computed  from  the  rifling  twist 
as  defined  in  the  stability  analysis  and  the  muzzle  velocity.  If  the  muzzle 
velocity  is  changed,  a  new  spin  rate  will  be  computed. 

Lines  four  through  twelve  define  the  initial  atmospheric  conditions 
and  gun  projectile  exit  conditions. 

Line  thirteen  is  the  integration  time  step.  For  a  2  DOF  analysis  the 
program  will  default  to  0.02  second.  For  a  6  DOF  analysis  the  program  will 
use  the  value  stored  in  the  projectile  data  file.  For  an  initial  analysis  a 
value  of  zero  should  be  entered.  This  will  cause  the  program  to  compute  a 
value  equal  to  l/20th  the  nutation  frequency.  This  will  cause  an  automatic 
redisplay  of  the  menu.  The  time  step  may  now  be  modified  to  a  round  value  if 
desired. 


The  trajectory  analysis  will  run  until  one  of  three  conditions  is 
reached;  the  projectile  strikes  ground  (z<*0),  the  maximum  time  is  reached 
(line  14),  or  the  maximum  slant  range,  the  distance  from  the  muzzle  to  the 
projectile  in  3D  space,  is  reached  (line  15).  If  the  projectile  strikes  the 
ground,  the  program  will  interpolate  the  output  data  back  to  z*0. 

Trajectory  output  may  be  created  by  one  of  two  methods,  even 
increments  of  time  or  slant  range.  The  program  will  allow  a  maximum  of  500 
output  points  to  be  calculated.  Caution  should  be  taken  in  using  slant  range 
output  with  high  angles  of  attack.  As  a  projectile  comes  over  the  top  of  the 
trajectory  it  is  possible  for  the  slant  range  to  start  decreasing. 

If  the  user  selects  to  examine/modify  the  aerodynamic  coefficients,  Table  8 
will  be  displayed.  The  look  of  this  menu  will  depend  on  the  type  of 
projectile  and  terminal  being  used.  For  a  fin  stabilized  projectile  there  are 
a  total  of  three  menus.  A  spin  stabilized  projectile  has  two.  In  the  small 
screen  mode  these  menus  each  take  two  screens  to  display.  To  make 
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TABLE  8.  AERODYNAMIC  COEFFICIENTS 
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oy  i»  "DONS*  t a  continue 
oy  in  “LIST"  to  rolist" 

oy  in  • Fillt'  for  i  Lino  prlntor  output 


notifications  to  these  tables  first  enter  the  column  number.  Next,  enter  line 
numbers  and  new  values  for  each  Mach  number  to  be  modified.  Enter  a  <er>  with 
no  line  number  when  modifying  a  column.  If  an  "S"  is  entered  after  the  <cr>, 
all  values  in  the  column  with  Mach  numbers  higher  than  the  last  one  modified 
will  be  given  the  values  last  entered. 

If  a  fumer  analysis  was  selected  on  menu  number  1,  the  following  Fumer 
Input  Menu  will  be  displayed: 

Trajectory  -  FUMER  Analysis  -  Input 


1  -  Axial  Inertia  at  Burnout  (lb*-in**2): 

2  -  Transverse  Inertia  at  Burnout  ( lbm-in**2 ) : 

3  -  Projectile  Height  at  Burnout  (lba): 

4  -  Center  of  Gravity  (calibers  fron  nose): 

5  -  Multiple  No sale  Cant  Angle  (deg): 

6  -  Multiple  Nozzle  Exit  Radius  (in): 

(fron  projectile  centerline) 

7  -  Burn  Start  Tine  (sec): 

8  -  Burn  Stop  Tine  (sec): 

9  -  Do  you  «rish  to  input  the  drag  vs  tine  : 

profile  or  the  drag  vs  Mach  no.  curve? 

(0:Continue  l:Drag  &  Thrust  vs  Tine  2:Drag  vs  Mach  No.) 

Key  in  iten  nunber  new  value 

Key  in  "LIST"  to  relist  "DONE"  to  continue 

Key  in  "PRINT"  for  a  Line  printer  output 

On  this  menu  the  physical  properties  of  the  projectile  after  fumer 
burn  out  are  defined.  The  Physical  Properties  segment  computes  properties 
with  and  without  the  tracer  elements.  Properties  are  assumed  to  vary  linearly 
during  the  time  of  the  tracer/fumer  burn.  It  is  not  necessary  that  tracer 
elements  exist  to  perform  a  fumer  analysis.  If  more  than  one  nozzle  exists, 
the  radius  of  the  nozzle  cluster  and  cant  angle  of  the  nozzles  may  be 
specified.  The  program  assumes  the  nozzles  are  arranged  in  a  circular 
pattern.  The  values  entered  for  thrust  are  the  thrust  along  the  nozzle 
orientation.  The  thrust  along  the  projectile  axis  will  be  the  thrust  times 
the  cosine  of  the  cant  angle.  A  rotational  force  of  thrust  times  the  sine  of 
the  cant  angle  is  also  computed  for  each  nozzle. 

Line  nine  allows  for  the  specification  of  projectile  thrust,  and/or 
drag  reduction  during  the  projectile  burn.  An  entry  of  one  or  two  will  cause 
additional  input  menu(s)  to  be  displayed.  An  entry  of  one  will  display: 

Trajectory  -  FUMER  Analysis  (Drag  &  Thrust  vs  Time)  -  Input 


0.007581 

0.058568 

0.230478 

3.560222 

0.000000 


0.000000 

2.400000 

2 


Decrease  in  Drag  Coef.  as  a  Function  of  Time 
Entry  #  Burn  Time  Delta  CX  Thrust 

sec  lbf 

Key  in  "DONE"  to  continue 
Key  in  "LIST"  to  relist 
Key  in  "PRINT"  for  a  Line  printer  output 
Key  in  values  separated  by  commas  to  modify 
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Up  to  20  lines  of  input  may  be  specified  here.  Each  line  must  contain 
a  time,  delta  drag,  and  thrust  entry.  Zero  may  be  entered  for  thrust  and 
drag.  Time  must  be  entered  in  ascending  order.  The  drag  reduction  entered 
here  is  magnitude  reduction  (negative  for  increase).  During  the  analysis,  the 
program  will  perform  a  linear  interpolation  from  this  table. 

Option  Two  in  line  nine  allows  for  £he  specification  of  a  reduced  drag 
versus  Mach  number.  A  table  of  drag  versus  Mach  number  with  and  without  fumer 
burn  will  be  displayed.  If  stability  analysis  has  been  run,  and  this  is  a 
spin  stabilized  projectile,  a  base  drag  from  McDrag  will  be  displayed. 
Otherwise  the  base  drag  is  specified  the  same  as  the  total  drag.  The  base 
drag  is  displayed  only  for  reference  purposes.  The  analysis  only  uses  the 
total  drag: 


Trajectory  -  FUMER  Analysis  (Drag  vs  Mach  Number)  -  Input 


No  Funer  During  Burn 

Mach  Humber  Total  Drag  Base  Drag  Total  Drag  Base  Drag 


0.0100 

0.1980 

0.1384 

1: 

0.1822 

18: 

0.1384 

0.6000 

0. 1980 

0.1411 

2: 

0.1822 

19: 

0.1411 

0.8000 

0.2004 

0.1477 

3: 

0.1844 

20: 

0.1477 

0.9000 

0.2200 

0.1645 

4: 

0.2024 

21: 

0.1645 

0.9500 

0.2590 

0.1776 

5: 

0.2383 

22: 

0.1776 

1.0000 

0.3710 

0.2518 

6: 

0.3413 

23: 

0.2518 

1.0500 

0.4270 

0.2122 

7: 

0.3928 

24: 

0.2122 

1.1000 

0.4180 

0.2169 

8: 

0.3846 

25: 

0.2169 

1.2000 

0.3920 

0.2108 

9: 

0.3606 

26: 

0.2108 

1.3500 

0.3670 

0.2015 

10: 

0.3376 

27: 

0.2015 

1.5000 

0.3460 

0.1915 

11: 

0.3183 

28: 

0.1915 

1.7500 

0.3130 

0.1733 

12: 

0.2880 

29: 

0.1733 

2.0000 

0.2820 

0.1543 

13: 

0.2594 

30: 

0.1543 

2.5000 

0.2300 

0.1200 

14: 

0.2116 

31: 

0.1200 

3.0000 

0.1900 

0.0927 

15: 

0. 1748 

32: 

0.0927 

4.0000 

0.1520 

0.0631 

16: 

0.1398 

33: 

0.0631 

5.0000 

0. 1300 

0.0454 

17: 

0.1196 

34: 

0.0454 

35  Total  burn  drag  as  Z  of  total  drag 

36  Base  burn  drag  as  Z  of  base  drag 

Key  in  item  number  new  value 

Key  in  "LIST”  to  relist  "DONE"  to  continue 

Mote  that  modifications  will  be  reflected  against  each  other 

Key  in  "PRINT"  for  a  Line  printer  output 

This  table  includes  four  options  for  modifying  the  drag.  The  total 
drag  or  base  drag  during  burn  may  be  changed  directly.  This  is  done  by 
modifying  lines  one  through  thirty-four. 
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For  any  change  to  total  drag,  the  same  magnitude  change  to  base  drag 
at  that  Mach  number  will  be  computed.  Conversely,  changes  in  base  drag  are 
reflected  on  the  total  drag.  Options  Thirty  Five  and  Thirty  Six  allow  for 
across  the  board  changes  by  a  specified  percentage.  With  Option  Thirty  Five, 
a  fractional  value  is  specified  by  which  the  total  drag  for  all  Mach  numbers 
will  be  multiplied.  The  magnitude  change  will  then  be  subtracted  from  the  base 
drag.  Option  Thirty  Six  performs  the  same  operation  on  the  base  drag.  The 
fractional  value  used  is  stored  in  the  projectile  model  file  and  will  be 
displayed  when  this  operation  is  selected.  It  is  displayed  only  for  reference 
purposes  so  it  must  be  reentered. 

If  Options  Thirty  Five  or  Thirty  Six  are  selected  the  table  will  be 
redisplayed  in  the  modified  form  automatically. 

If  desired,  both  the  drag  versus  time  and  drag  versus  Mach  number 
tables  may  be  used.  For  the  case  of  a  rocket  it  will  be  necessary  to  use  the 
drag  versus  time  menu  to  specify  the  projectile  thrust  profile.  Depending  on 
how  the  drag  reduction  data  is  specified  it  may  be  necessary  to  use  either  or 
both  menus  to  define  the  drag  reduction. 

If  an  aircraft-originated  trajectory  was  specified  in  menu  number  one, 
the  following  menu  will  be  displayed  for  specification  of  the  aircraft 
parameters: 


Trajectory  -  Aircraft  Originated  —  Input 


Aircraft  Inputs 

1  -  Aircraft  Velocity  (ft  /  sec):  0.000 

2  —  Aircraft  Dive  Angle  (deg.,  +”down):  0.000 

3  -  Aircraft  Dive  Azimuth  (deg.,  +“nose  left):  0.000 

4  —  Aircraft  Pitch  Angle  (deg.,  +*nose  down):  0.000 

5  -  Aircraft  Yaw  Angle  (deg.,  +“nose  left):  0.000 

6  —  Gun  Elevation  WRT  Aircraft  (deg.,  +“up):  0.000 

7  -  Gun  Azinuth  WRT  Aircraft  (deg.,  +*left):  0.000 


Key  in  iten  number  n,N  new  value 

Key  in  "LIST"  to  relist  "DONE"  to  continue 

Key  in  "PRINT"  for  a  Line  printer  output 

For  an  aircraft-originated  trajectory,  the  elevation  angle  used  for 
the  analysis  will  be  the  dive  angle  from  this  menu.  Any  other  value  entered 
will  be  ignored. 
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If  a  ball  rotor  simulation  was  selected  on  menu  number  one,  the 
following  menu  will  be  displayed: 

Trajectory  -  Ball  Rotor  Simulation  -  Input 

*  -  -  -  -  - - - - -  -  -  - 


1  -  Ball  Height  (grans):  0.0000 

2  —  Ball  CG  from  Projectile  CG  (in):  0.0000 

3  —  Ball  Clearance  (in):  0.0000 

Key  in  item  number  ","  new  value 

Key  in  "LIST"  to  relist  "DONE"  to  continue 

Key  in  "PRINT"  for  a  Line  printer  output 

If  more  than  one  barrel  was  specified  in  menu  number  one,  the  Multiple 
Barrel  Input  menu  will  be  displayed: 

Trajectory  -  Multiple  Barrel  Gun  —  Input 


Multibarrel  Gun  Inputs 


Gun  Half— cone  Angle 

(deg): 

0.000 

Muzzle  Radius 

(in): 

0.000 

Sear-off  Position 

(deg,  +*cdu): 

0.000 

from  bottom  vertical  looking  downrange 

Rate  of  Fire 

(SPM): 

0.000 

Action  Time 

(msec): 

0.000 

Key  in  item  number  new  value 

Key  in  "LIST"  to  relist  "DONE"  to  continue 

Key  in  "PRINT"  for  a  Line  printer  output 


These  inputs  define  the  rotational  motion  of  a  Gatling  gun  barrel 
cluster.  This  data  will  determine  the  angular  forces  applied  to  the 
projectile  on  exit  from  the  barrel. 


This  completes  the  data  input  for  trajectory  analysis.  The  program 
will  perform  the  analysis  and  integrate  the  equations  of  motion.  A  fourth 
order  Runge-Kutta  integration  is  used.  The  results  of  the  last  integration 
step  will  be  printed  out  followed  by  the  output  options: 


Terminal  Parameters 


Velocity  *  690.244  (ft/sec) 

X-dist  -  33119.441  (ft) 

Y-dist  «  -129.914  (ft) 

Z-dist  -  0.000  (ft) 

Time  -  24.48492  (sec) 
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Once  the  terminal  parameters  have  been  displayed,  the  output  option 
menu  will  appear: 

6  DOF  Trajectory  Segment  Output 
Plots 

Time  of  Flight(sec)  Slant  Range(ft) 


1 

X  (ft) 

21 

2 

Y  (ft) 

22 

3 

Z  (ft) 

23 

4 

Slant  Range  (ft) 

24 

5 

Velocity  (ft/sec) 

25 

6 

Mach  Number 

26 

7 

Spin  Rate  (rad/sec) 

27 

8 

Total  Yaw  (deg) 

28 

9 

Drop 

29 

10 

Drift 

30 

11 

Gyro.  Stability  Factor 

31 

12 

Pitch  [THETAj-M  (deg) 

32 

13 

Yaw  [PS1]-M  (deg) 

33 

14 

Angular  Accel  (rad/sec2) 

34 

15 

Angular  Rate  (rad/sec) 

35 

16 

Dynamic  Stability  Factor 

36 

17 

Kinetic  Energy 

37 

38  Y  vs  X 

39  Z  vs  X 

40  Theta  vs  Psi 

50  Return  to  Input 

51  Cross  Plots 

52  Return  to  Main  Menu 

54  Tic  nark  switch  (on) 

Tables 

61  Complete  Output 

62  T  vs  X,Y,Z,XBAR,V, Spin, Drop, Deflection 

63  T  vs  X,Y,Z, Phi, Theta, Psi, Alpha 

64  T  vs  X,Y,Z  for  Projectile  WRT  Barrel  C.L. 

Gun  C.L. ,  Aircraft  C.L. 

65  T  vs  X,Y,Z,V,RE 

Plots  are  against  either  time-of-f light  or  slant  range.  Tables  are 
computed  for  time-of-f light  and  will  be  sent  to  the  designated  output  printer. 

Gyroscopic  and  dynamic  stability  plots  are  only  available  if  the 
linear  theory  option  was  selected  (input  menu  number  two). 

Pitch  and  yaw  data  is  only  available  if  a  6  DOF  analysis  was  run. 
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Option  Fifty  One  will  send  the  program  to  the  cross  plot  menu.  The 
names  and  descriptions  of  the  trajectories  available  for  plotting  will  be 
listed.  Up  to  10  trajectories  may  be  selected  from  this  list  for  plotting 
together.  Approximately  500  points  will  be  plotted  together.  If  each  curve 
contains  500  points  and  three  curves  are  to  be  plotted  together,  then  every 
third  point  will  be  plotted.  The  more  curves  that  are  plotted  together,  the 
fewer  points  on  each  curve. 

Option  Fifty  Four  will  toggle  the  tick  mark  flag.  The  option  specifies 
which  way  the  toggle  will  go.  For  a  slant  range  plot  the  tick  marks  will  show 
each  second  of  time  of  flight.  For  a  slant  range  plot,  the  tick  marks  will 
show  each  1000  feet/meters  of  slant  range. 

6  DOF  Trajectory  Cross  Plot  Segment 


Time  of  Flight  Slant  Range(ft) 

1  X  (ft)  21 

2  Y  (ft)  22 

3  Z  (ft)  23 

4  Slant  Range  (ft)  24 

5  Velocity  (fps)  25 

6  Mach  No.  26 

7  Spin  late  (rad/sec)  27 

8  Total  Yaw  (deg)  28 

9  Drop  (arad)  29 

10  Drift  (arad)  30 

11  Gyroscopic  Stability  Factor  31 

12  Theta*-M  (deg)  32 

13  Psi-M  (deg)  33 

14  Dynaaic  Stability  Factor  34 

15  Kinetic  Energy  35 


38  Y  vs  X 

39  Z  vs  X 

50  Return  to  Input 

52  Return  to  Main  Menu 

53  Return  to  Plot  Menu 

54  Tic  aark  switch  (on) 

Figure  16  shows  a  cross  plot  of  Z  vs  X  for  elevation  angle  variation. 
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Figure  16.  Cross  Plotted  Trajectory 


SECTION  VI 


LINEAR  THEORY 


1.  GENERAL 

Within  one  executable  program  both  the  data  input  and  analysis  are 
performed.  Once  data  has  been  input  it  may  be  stored  for  future  reference. 
The  following  options  menu  will  be  displayed  for  data  input/analysis: 

Enter  the  number  of  the  desired  option 

1  -  Edit  Input  Data 

2  -  Run  Analysis 

3  -  Store  Master  File  Data 

4  —  Exit  Linear  Theory  Program 


The  program  will  not  allow  for  analysis  (Option  Two)  until  data  has  been 
returned  from  the  editing  module  (Option  One).  Data  may  be  stored  (Option 
Three)  before  or  after  the  running  of  an  analysis.  It  is  advised  that  data  be 
stored  on  exiting  the  edit  module.  Storing  input  at  this  time  will  prevent 
data  being  lost  if  a  hard  crash  occurs  during  the  analysis.  If  the  results  of 
the  analysis  are  acceptable,  it  is  advised  that  the  data  be  stored  again.  By 
doing  this,  the  parameters  found  from  the  analysis  will  be  the  initial 
estimates  for  any  future  reduction  of  this  shot. 

2.  LINEAR  THEORY  INPUT 

Upon  entering  the  data  Edit/Input  module,  the  following  menu  of  options 
will  be  displayed: 

1  -  Initialize  new  data  set 

2  -  Retrieve  cataloged  data  set 

3  -  Shot  number  and  description 

4  -  Linear  Theory  analysis  Flags 

5  -  Projectile  physical  properties 

6  —  Analysis  reference  data 

7  -  Initial  aerodynamic  coefficient  estimates 

8  -  Initial  aerodynamic  parameters  estimates 

9  -  Aerodynamic  fit  flags 

10  --  Atmospheric  parameters 

11  -  Edit  current  Tunnel  XYZ  data 

12  -  Use  new  data  set  from  Tunnel  XYZ  file 

13  -  Delete  current  Tunnel  XYZ  data 

14  -  Done 
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Each  of  the  above  options  will  cause  a  different  menu  of  parameters  to  be 
presented  for  editing.  A  common  procedure  for  editing  has  been  established. 
The  line  number  and  new  value  are  entered*  with  one  change  per  line.  A  "LIST" 
option  will  redisplay  the  data  with  the  changes.  "DONE"  will  terminate  the 
edit  session  and  return  to  the  above  menu.  "R"  will  cause  the  program  to  step 
through  each  entry  in  the  edit  menu.  The  line  numbers  within  the  menu  will  be 
incremented  automatically. 

The  first  step  in  an  editing  session  is  either  to  retrieve  an  existing 
data  set  or  to  initialize  a  new  data  set.  Edit  Option  One  is  used  to 
initialize  a  new  data  set.  Under  this  option*  the  program  will  automatically 
step  through  Options  Three  through  Ten  and  Option  Twelve. 

To  retrieve  data  stored  during  a  previous  data  reduction  session  use  Edit 
Option  Two.  If  the  user  chooses  to  retrieve  existing  data,  a  question  will 
appear  requesting  the  shot  number.  A  list  of  available  shots  may  be  displayed 
at  this  time: 


Key  in  10  digit  projectile  title  (or  LIST):  LIST 

Enter  the  shot  group  number  to  list:  -1 

Linear  Theory  Projectile  Titles 
Previous  Reduction 

Shot 

Group  Shot  Number  Date  Time  Description 


9 

9 

9 

8 

8 

0 

0 

0 

9 


BS87040819 

BS87040613 

BS87040209 

BS87080617 

BS87080618 

BS87032490 

BS87032797 

BS87042434 

BS87042333 


1-SEP— 87 
1 -SEP-8 7 
l-SEP-87 
24— AUG-87 
24-AUG-87 
26— JUN-87 
26— JUN-87 
15-  JUN-87 
15— JUN-87 


13:04:00 

13:04:00 

11:57:20 

15:24:13 

15:24:13 

10:37:12 

10:37:12 

10:49:38 

10:49:38 


40  MM  HEDP  TUBULAR  REREAD 
40  MM  HEDP  TUBULAR  REREAD 
40  MM  HEDP  TUBULAR  REREAD 
20  MM  BOOMED /CASED 
20  MM  BOOMED/CASED 
40  MM  HEDP  TUBULAR 
40  MM  HEDP  TUBULAR 
20  MM  BOOMED/CASED 
20  MM  BOOMED/CASED 


Key  in  10  digit  projectile  title  (or  LIST):  BS87040716 


Note  that  in  the  above  table  a  negative  shot  group  number  was  entered.  A 
positive  group  number  will  cause  only  that  shot  group  to  be  displayed.  A 
negative  group  number  will  cause  all  shot  groups  to  be  displayed.  Within  this 
table  are  the  date  and  time  that  the  shot  was  previously  cataloged. 

Edit  Option  Three  will  allow  the  user  to  specify  a  new  shot  number  and 
description. 
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Edit  Option  Four  contains  the  analysis  option  flags: 

Linear  Theory  Analysis  Flags 


1 

2 

3 

4 

5 

6 


7 

8 


9 


10 


Starting  point  to  fit 

Length  of  first  section  to  be  fit 

Increaental  section  increment 

Stoning  section  increment 

Frequency  guess 

0:  Do  not  use  Caa  1:  Use  Caa 

Foil  Flag 

0:  Roll  data  not  available  (or  do  not  use) 
1:  Roll  data  available  (autoaatic  unwind) 
2:  Use  roll  data  as  is  (already  unwound) 

3:  Unwind  roll  data 
Autoaatic  guess  routine  flag 
0:  Do  not  use  1:  Use 

Spin  Flag 

0:  Spin  stabilized  (fixed  plane) 

1:  Fin  stabilized  (fixed  plane) 

2:  Fin  stabilized  (body  fixed) 

3:  Fin  stabilized  Caa  O  -CNB  (body  fixed) 

Sections  Flag 

0:  Use  LN,  NA,  ND,  NB 


ILH]:  1 

[NA] :  25 

[ND]:  0 

[NB] :  2 

[NAUTO]:  0 

[NR0LL] :  0 


[NSG] :  0 

[NSP1N ] :  0 


[NSECT]:  3 


1:  Use  autoaatic  sections  (Output  sectional  fits) 

2:  Use  autoaatic  sections  (Output  1  to  NT) 

3:  Help 

Shot  Group  Nuaber  [IGRP]:  0 


Key  in  itea  nuaber  new  value 

Key  in  "LIST"  to  relist  "DONE"  to  continue 

Key  in  "R"  to  step  through  input 


Within  this  table,  LN,NA,ND,NB  are  used  to  define  the  section  of  range 
stations  to  be  used  in  fitting  the  data  and  how  the  data  will  be  incremented. 


NR0LL  is  used  by  the  program  to  flag  the  routines  used  to  analyze  the  roll 
data. 


If  NSG  is  set  on  (set  to  1),  the  program  will  automatically  select  the 
parameters  to  fit. 

NSPIN  specifies  the  type  of  projectile,  fin. or  spin  stabilized,  and  the 
coordinate  system,  fixed  plane  or  body  fixed. 

If  NSECT  is  set  to  three,  the  HELP  routines  will  be  called.  These 
routines  give  the  user  a  chance  to  control  fit  curve  fitting  interactively  to 
help  get  things  started.  During  the  HELP  phase  only  yaw  motion  is  being 
analyzed.  Plotted  and  tabulated  data  is  available  for  examination.  (See  next 
section  on  Linear  Theory  Help.) 
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The  shot  group  number  is  initially  set  by  the  TDBIAS  program  and  may  be 
modified  here  if  desired. 

Edit  Option  Five  contains  the  physical  properties  of  the  projectile  to  be 
analyzed: 

Projectile  Properties 


1 

-  Projectile  diameter 

[D] 

(in): 

1.558 

2 

-  Axial  moment  of  inertia 

[IX] 

(slog-ft-in): 

0.0005250 

3 

-  T  axis  moment  of  inertia  [IT] 

(slug— ft-in): 

0. 

0010540 

4 

-  Z  axis  moment  of  inertia  [IZ] 

(slug— f t-in): 

0.0010540 

5 

-  Product  of  inertia 

[IXY] 

(slug-ft-in): 

0. 

0000000 

6 

-  Projectile  weight 

[WGT] 

(gram): 

244.77800 

7 

-  Center  of  gravity 

[ICG] 

(from  nose): 

0.612 

8 

-  Projectile  length 

[XL] 

(in): 

3.084 

Key  in  item  number  ","  new  value 

Key  in  "LIST"  to  relist  "DONE"  to 

continue 

Key  in  ' 

Tl"  to  step  through  input 

Option  Six  contains  the  analysis 

reference  parameters: 

Analysis  Reference  Parameters 

1  ~  Initial  roll  rate 

[CA] 

(deg/ft): 

93.00000 

2  -  Roll  pin  orientation 

[RBAIS ] 

(deg): 

0.00000 

3  -  Reference  mach  number 

[HREF] 

0.60000 

4  -  Reference  diameter 

[DIAR] 

(in)’: 

1.55764 

5  -  Reference  CG  location 

[XCGR] 

(from  nose): 

0.61200 

6  -  Reference  length 

[XLR] 

(in): 

3.08440 

7  -  Number  of  fins 

[XFIN] 

(symmetry  planes): 

1 

8  —  Roll  acceleration 

[PH2] 

(deg/ft**2): 

0.0000E+00 

Key  in  item  number  new  value 

Key  in  "LIST"  to  relist  "DONE"  to  continue 

Key  in  "R"  to  step  through  input 
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EdiC  Option  Seven  allows  for  modification  of  the  aerodynamic  coefficient 
estimates: 


Estimates  of  Aerodynamic  Coefficients 


rrrrt 


2  -  Cma  1.40000 

3  -  Clp  :  -0.02600 

4  —  Cld  :  0.00000 

5  -  CX  Mach  :  0.00000 

6  -  CXa2  :  2.30000 

7  -  Gan  Barrel  Twist  (cal/rev):  29.71000 


Key  in  item  number  new  value 

Key  in  "LIST"  to  relist  "DONE"  to  continue 
Key.  in  "1"  to  step  through  input 

Within  linear  theory  only  the  gun  barrel  twist  is  used.  The  twist  is  used 
to  determine  the  initial  roll  rate  of  the  projectile.  The  other  coefficients 
are  just  passed  through  the  data  file  for  initial  estimates  by  the  6  DOF 
program. 

Edit  Option  Eight  allows  for  modification  of  the  initial  aerodynamic 
parameter  estimates  used  by  linear  theory.  The  values  of  these  parameters 
will  change  during  the  linear  theory  reduction.  When  a  good  fit  is  achieved 
the  projectile  file  should  be  stored.  This  will  cause  the  final  (fit)  values 
of  the  estimates  to  be  retained  as  a  future  starting  point: 

Initial  Aerodynamic  Parameter  Estimates 


1  -  Nutation  vector 

[Kl] 

(deg): 

0.17710 

2  -  Precession  vector 

[K2] 

(deg): 

0.09820 

3  -  Nutation  damping  exponent 

IL1] 

(1/ft): 

0.00070 

4  -  Precession  damping  exponent 

[L2] 

(1/ft): 

0.00700 

5  ~  Nutation  vector  orientation 

[PI] 

(deg): 

299.06381 

6  -  Precession  vector  orientation 

[P2] 

(deg): 

273.14471 

7  -  Nutation  frequency 

[Wl] 

(1/ft): 

56.48310 

8  -  Precession  frequency 

[W2] 

(1/ft): 

23.84210 

9  -  Nutation  freq.  change 

(WD1  ] 

(deg/ft**2): 

0.00600 

10  -  Precession  freq.  change 

[WD2] 

(deg/ft**2): 

0.00340 

11  -  Trim  vector 

[K3J 

(deg): 

0.00000 

12  -  Trim  vector  orientation 

[P3] 

(deg): 

0.00000 

Key  in  item  number  new  value 

Key  in  "LIST"  to  relist  "DONE"  to  continue 

Key  in  "R"  to  step  through  input 
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Edit  Option  Nine  allows  for  modification  of  the  parameter  fit  flags.  The 
flags  that  are  set  determine  which  estimates  may  be  modified  by  the  program 
during  the  data  fitting: 

Aerodynamic  Modification  Fit  Flags 


1  -  Natation  vector  [Kl] 

2  -  Precession  vector  [K2] 

3  -  Natation  damping  exponent  [LI] 

4  -  Precession  damping  exponent  [L2] 

5  -  Natation  vector  orientation  [PI] 

6  -  Precession  vector  orientation  [P2] 

7  -  Natation  frequency  [HI] 

8  -  Precession  frequency  [W2] 

9  -  Nutation  freq.  change  [WD1] 

10  -  Precession  freq.  change  [WD2] 

11  -  Trim  vector  [K3] 

12  -  Trim  vector  orientation  [P3J 

(0  -  do  not  fit  :  1  -  fit) 

13  -  Restart  flag  [RES] 

(0  -  Restart  :  1  -  No  restart) 


Key  in  item  number  MSM  new  value 

Key  in  "LIST"  to  relist  "DONE"  to  continue 

Key  in  "R"  to  step  through  input 

If  the  HELP  option  is  to  be  run,  these  flags  may  be  changed  within  that 
module. 

Edit  Option  Ten  allows  for  modification  of  the  atmospheric  parameters. 
Either  relative  humidity  or  air  density  may  be  specified.  If  either  value  is 
changed,  the  other  will  be  computed  and  displayed: 

Atmospheric  Conditions 


1  ~  Air  temperature  (Deg  C): 

2  -  Atmospheric  pressure  (mbar): 

3  -  Relative  humidity  : 

or 

4  -  Air  density  (slug/ft**3): 

Key  in  item  number  new  value 
Key  in  "LIST"  to  relist  "DONE"  to  continue 


19.600 

1026.200 

0.460 

0.002361 
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Edit  Option  Eleven  allows  the  user  to  modify  time  and  roll  data  for  any  of 
the  stations  along  the  range.  A  table  of  the  time  and  roll  data  for  each 
rrange  station  along  with  the  station  distances  will  be  displayed.  Station 
distances  may  not  be  changed,  they  are  only  displayed  for  reference  purposes. 
Entering  a  negative  time  will  remove  a  station  from  the  analysis: 

Tunnel  ZXZ  data  for  editing 


Enter  negative  time  to  remove  a  station 


Station 

Distance 

(feet) 

Time 

(sec) 

Roll 

(deg) 

Station 

Distance 

(feet) 

Time 

(sec) 

Roll 

(deg) 

11.9377 

(  1) 

0.5584267 

(14) 

-99. 

242.0256 

(27) 

0.7894410 

(39) 

-99. 

17.1983 

(  2) 

0.5634965 

(15) 

-99. 

252.0643 

(28) 

0.7999383 

(40) 

-99. 

27.1724 

(  3) 

0.5731447 

(16) 

-99. 

271.9348 

(29) 

0.8207810 

(41) 

-99. 

47.2351 

(  4) 

0.5925861 

(17) 

-99. 

276.9530 

(30) 

0.8260718 

(42) 

-99. 

57.1083 

(  5) 

0.6022863 

(18) 

-99. 

291.9944 

(31) 

0.8419708 

(43) 

-99. 

67.0824 

(  6) 

0.6120735 

(19) 

-99. 

306.9884 

(32) 

0.8578804 

(44) 

-99. 

72.1105 

(  7) 

0.6170182 

(20) 

-99. 

322.0588 

(33) 

0.8739316 

(45) 

-99. 

76.9907 

(  «) 

0.6218238 

(21) 

-99. 

352.1425 

(34) 

0.9062108 

(46) 

-99. 

151.9478 

(  9) 

0.6967998 

(22) 

-99. 

382.1923 

(35) 

0.9387400 

(47) 

-99. 

181.9187 

(10) 

0.7273291 

(23) 

-99. 

397.1584 

(36) 

0.9550332 

(48) 

-99. 

186.9524 

(ID 

0.7324781 

(24) 

-99. 

412.1839 

(37) 

0.9714601 

(49) 

-99. 

212.0999 

222.1359 

(12) 

(13) 

0.7583820 

0.7687663 

(25) 

(26) 

-99. 

-99. 

442.1961 

(38) 

1.0044963 

(50) 

-99. 

Key  in  item  number  new  value 

Key  in  "LIST"  to  relist  "DONE"  to  continue 

If  a  station  has  been  removed  a  -99  will  be  displayed  for  the  roll  angle. 

Edit  Option  Twelve  allows  for  the  replacement  of  the  Tunnel  XYZ  data  set 
with  the  data  for  another  shot  number.  The  projectile  name  and  description 
will  be  replaced  as  well,  with  that  of  the  new  Tunnel  XYZ  data  set.  All  other 
data  will  remain  the  same  as  it  was. 

Edit  Option  Thirteen  removes  all  Tunnel  XYZ  data  but  all  other  data 
remains . 

3.  LINEAR  THEORY  HELP 

If  the  Linear  Theory  HELP  option  is  selected  in  the  Linear  Theory  Input 
(setting  flag  number  nine  to  three  in  Edit  Option  Four),  the  program  will 
branch  into  the  HELP  mode  and  remain  there  until  directed  otherwise.  This 


» 
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nodule  performs  a  reduction  on  the  yaw  data  and  helps  to  refine  the  initial 
estimates  before  running  the  full  linear  theory  reduction.  The  HELP  input 
menu  is: 

BS87040716  40ms  HEDP  Tabular 
28-5EP— 87  16:18:26 

Station  Options 


1  ~  Starting  Point  to  Fit  :  10  11.94  (ft) 

2  -  Length  Of  First  Fit  Section  :  25  6  442.20  (ft) 

3  -  Summing  Section  Increment  :  2 


Initial  Estimates  Fit  Flag 


11 

— 

1.0000 

Nutation  Vector 

[Kl] 

(deg): 

1 

-  31 

12 

- 

1.0000 

Precession  Vector 

[K2J 

(deg): 

1 

-  32 

13 

- 

0.0007 

Nutation  Damping  Exponent 

[LI] 

(1/ft): 

0 

-  33 

14 

— 

0.0070 

Precession  Itemping  Exponent 

[L2] 

(1/ft): 

0 

-  34 

15 

- 

299.0000 

Natation  Vector  Orientation 

[PI] 

(deg): 

1 

-  35 

16 

- 

273.0000 

Precession  Vector  Orientation 

[P2] 

(deg): 

1 

-  36 

17 

- 

44.0000 

Nutation  Frequency 

[HU 

(1/ft): 

1 

-  37 

18 

- 

25.0000 

Precession  Frequency 

[H2] 

(1/ft): 

1 

-  38 

19 

— 

0.0060 

Natation  Frequency  Change 

[WD1] 

(deg/ft**2): 

0 

-  39 

20 

— 

0.0034 

Precession  Frequency  Change 

[WD2] 

(deg/ft**2): 

0 

-  40 

21 

- 

0.0000 

Trim 

[«] 

(deg): 

0 

-  41 

22 

— 

0.0000 

Trim  Orientation 

[P3] 

(deg): 

0 

-  42 

Key  in  item  number  new  value 

Key  in  "LIST"  to  relist  "DONE"  to  continue 


The  distances  shown  in  lines  one  and  two  are  station  locations  and  are 
output  for  reference  purposes.  Once  the  input  data  is  satisfactory,  enter 
DONE  to  run  the  analysis. 

The  HELP  output  will  display  the  solution  for  the  parameters  and  the 
probable  error  in  the  yaw: 

BS87040716  40mm  HEDP  Tubular 
28-SEP— 87  16:18:26 

HELP  Results 


Yaw  Error  (deg) 

Starting  Point  to  Fit 
Length  Of  First  Fit  Section 
Summing  Section  Increment 
Total  Number  of  Sections 


0.314 


1 

11.94 

(ft) 

25 

442.20 

(ft) 

2 

25 

442.20 

(ft) 
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BS87Q40716  40a  HEDP  Tubular  -  (Continued) 


Analysis  Parameters 


Nutation  Vector  [XI]  (deg):  0.39889 
Precession  Vector  [K2]  (deg):  I .051 30 
Natation  Daaping  Exponent  [LI]  (1/ft):  0.00070 
Precession  Damping  Exponent  [L2]  (1/ft):  -0.00308 
Nutation  Vector  Orientation  [PI]  (deg):  288.82666 
Precession  Vector  Orientation  [P2]  (deg):  274.06430 
Natation  Frequency  [VI]  (deg/ft):  44.45790 
Precession  Frequency  [W2]  (deg/ft):  25.36411 
Imitation  Frequency  Change  [WD1]  (deg/ft**2):  -0.00067 
Precession  Frequency  Change  [WD2]  (deg/f t**2):  0.02859 
Trim  [E3]  (deg):  0.00000 
Trim  Orientation  [P3]  (deg):  0.00000 


At  this  point,  one  of  the  output  options  may  be  selected.  A  <cr>  will 
cause  the  option  menu  to  be  displayed  if  the  number  desired  is  not  known. 

The  following  output  options  exist  within  the  HELP  module: 

Select  desired  option 

1  -  Add  section  increment  and  continue  in  Help 

2  -  Edit  input  parameters 

3  -  Reset  estimates  to  original 

4  -  Reset  estimates  to  previous  analysis 

5  -  Display  data  list  (1-MT) 

6  -  Display  yaw  fit  (LN-NA) 

7  -  Exit  Help  /  Resume  Linear  Theory  analysis 

8  —  Exit  analysis  program 

9  -  Plot  Alpha  vs  X 

10  —  Theta  vs  X 

11  —  Psi  vs  X 

12  -  Theta  vs  Psi 

13  -  Theta  &  PSI  vs  X 

14  -  Theta  vs  Psi  (step) 


Option  One  will  increment  the  length  of  the  fit  by  the  number  specified  in 
•  the  summing  section  increment.  The  HELP  routine  will  then  be  repeated  for 

this  set  of  conditions. 
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Option  Two  will  return  the  program  to  the  input  menu  in  HELP.  At  this 
point  the  user  may  change  parameter  estimates  or  fit  flags.  The  program  will 
then  analyze  the  motion  for  the  new  set  of  conditions. 

Option  Three  causes  the  fit  flags  and  parameter  estimates  to  be  returned 
Co  what  they  were  when  the  HELP  module  was  entered.  The  program  will  proceed 
to  the  input  edit  menu. 

Option  Four  will  reset  the  estimates  to  what  they  were  before  the  last 
pass  through  the  fitting  routines.  The  program  will  return  to  the  edit  menu 
for  user  action. 


Option  Five  will  display  a  table  of  the  motion  data  as  read  from  the  film. 
The  following  shows  a  partial  list  of  this  output: 

Down-Range  Horizontal  Vertical 


Time 

Travel 

[X] 

Motion 

IT] 

Motion 

[z] 

Pitch 

[THETA] 

Taw  1 

[PSI] 

Roll  law  Roll 
[PHI] 

(sec) 

(ft) 

(feet) 

(feet) 

(deg) 

(deg) 

(deg) 

(deg) 

0.0000000 

11.9377 

4.56981 

0.91997 

-0.04756 

-1.02851 

0.0 

0.0 

0.0050698 

17.1983 

4.59496 

0.99552 

-0.15527 

0.51796 

0.0 

0.0 

0.0147180 

27.1724 

4.64285 

1.10731 

-0.23319 

-0.29336 

0.0 

0.0 

0.0341594 

47.2351 

4.80058 

1.32640 

-0.91906 

5.79866 

0.0 

0.0 

0.0438596 

57.1083 

4.79740 

1.43508 

-0.60334 

0.27905 

0.0 

0.0 

0.0536468 

67.0824 

4.85454 

1.54357 

-0.93052 

-1.11118 

0.0 

0.0 

Option  Six 

displays 

a  list  of 

the  yaw 

fit  together  with 

the  pitch  and 

from  the  film  reading.  The  following  shows  a  partial  display  from  a 
reduction: 


Time  Travel  [X] 
(sec)  (feet) 


0.0050698 

0.0147180 

0.0341594 

0.0438596 

0.0536468 


11.938 

17.198 

27.172 

47.235 

57.108 

67.082 


Pitch-Missile  Taw-Missile 

Experimental  Calculated  Experimental  Calculated 


(deg) 

(deg) 

(deg) 

(deg) 

0.775 

1.347 

-1.302 

-1.173 

0.549 

0.788 

0.243 

-0.160 

0.398 

0.527 

-0.686 

-0.251 

-0.292 

-0.106 

5.503 

9.754 

0.024 

-1.123 

0.123 

0.423 

-0.311 

-0.656 

-1.445 

-2.412 

Option  Seven  will  terminate  the  HELP  function  and  continue  Linear  Theory 
with  the  current  set  of  parameter  estimates. 

Option  Eight  will  terminate  the  Linear  Theory  analysis. 


Options  Nine  through  Fourteen  will  create  plots  of  the  projectile  motion 
as  determined  by  the  last  pass  through  the  data  reduction.  The  experimental 
points  will  be  superimposed  on  the  motion  plots  for  the  length  of  the  fit 
being  analyzed.  Figure  17  shows  pitch  and  yaw  data  for  travel  down  the  range. 
Figure  18  shows  a  Theta  versus  Psl  step  plot.  On  the  step  plots  the  points 
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BS87040716  40  MM  HEDP  TUBULAR 
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Figure  17.  Linear  Theory  Angular  Motion  Fit 


LINEAR  THEORY  REDUCTION: 


O 


Q.HI-UI  II—  X  UJ  I—  <C  I  O  LU  O 
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Figure  18.  Linear  Theory  Vector  Orientation 


are  displayed  one-at-a-time  with  a  pause  for  Che  user  to  type  <space  bar> 
between  points.  Some  terminals  aay  also  require  a  <cr>.  If  any  letter  is 
typed  other  than  <space  bar>  the  plot  will  terminate  there.  For  each  point 
displayed  in  the  step  plot  the  nutation  and  precession  vectors  are  drawn. 

4.  LINEAR  THEORY  ANALYSIS 

Linear  Theory  Analysis  aay  be  run  either  directly  froo  the  Input  module  or 
after  running  the  HELP  module.  The  user  will  be  prompted  for  a  number  of 
output  options. 

If  roll  data  exists  and  the  user  selected  the  unwind  option,  the  user  has 
the  option  of  examining  this  data. 

Once  the  analysis  has  been  performed,  the  user  may  examine  the  analysis 
summary.  It  is  recommended  that  this  data  be  exaained  before  .cataloging  the 
data. 

The  user  must  select  whether  or  not  to  catalog  the  data  in  the  sumnary 
file.  If  data  already  exists  for  this  shot  number,  it  will  be  overwritten. 

The  user  may  choose  to  use  the  same  output  results  for  plotting.  Plot 
data  is  maintained  by  shot  number  so  the  output  data,  once  saved,  will  be 
available  until  the  user  decides  to  delete  the  data. 

Two  forms  of  tabulated  output  are  available.  The  full  printout  and  a 
summary  printout.  The  summary  printout  is  a  subset  of  the  full  printout  and 
is  not  offered  for  printing  if  the  full  printout  has  already  been  selected. 
Table  9  shows  an  abbreviated  summary  printout. 

5.  LINEAR  THEORY  OUTPUT 

The  output  program  for  linear  theory  allows  the  user  to  create  plots  of 
the  coefficients,  experimental  points,  or  integrated  motion.  Tabulated 
summaries  may  be  made  of  any  shot  group. 

The  coefficient  plot  menu  is  as  follows: 

Summary  Segment  -  Linear  Theory  Reductions 


Shot  Group  0  -  40mm  HEDP 


♦ 


* 


1  -  CD 

2  -  CD  Mach 

3  -  CDO  Mach 

4  -  CHa 

5  -  Cma 

6  -  Cnpa 

7  -  Cmq 


99  -  Done 
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TABLE  9.  LINEAR  THEORY  ANALYSIS  SUMMARY 
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TABLE  9.  LINEAR  THEORY  ANALYSIS  SUMMARY  (CONTINUED) 
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TABLE  9.  LINEAR  THEORY  ANALYSIS  SUMMARY  (CONCLUDED) 
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Figure  19  shows  a  plot  of  the  linear  theory  drag  coefficients  from  a 
reduction  of  40mm  projectiles. 

The  following  menu  is  available  for  experimental  point  plotting: 
Linear  Theory  Experimental  Point  Plots 


1  -  Time 

2  -  Horizontal  Motion  [T] 

3  -  Vertical  Motion  [Z] 

4  -  Pitch  [THETA 1 

5  -  Taw  [PSI1 

6  -  Roll  [PHI] 

7  -  Pitch/Taw  [THETA-PSI] 

8  -  Pitch/Taw  Step  Plot 

99  -  Done 


The  experimental  points  will  be  superimposed  on  the  motion  plots.  A  plot 
of  experimental  data  was  shown  in  Figure  3. 

The  following  menu  is  available  for  the  motion  plots: 

Enter  the  number  for  the  desired  plot 

1  ~  Angle  of  Attack  [ALPHA]  vs  Travel  [X] 

2  -  Pitch  [THETA]  vs  Travel  [X] 

3  -  Yaw  [PSI]  vs  Travel  [X] 

4  -  Pitch  [THETA]  vs  Taw  [Psi] 

5  -  Theta  &  Psi  vs  X 
99  -  Done 


Figure  17  in  the  Linear  Theory  HELP  paragraph  shows  a  plot  of  the  same 
form  as  will  be  created  here.  The  code  for  plotting  in  the  HELP  option  is  the 
same  as  was  used  for  these  output  plots. 

Summary  tables  are  created  by  shot  group  number.  For  the  linear  theory 
program  the  summary  tables  are  sorted  by  both  shot  number  and  Mach  number. 
Table  10  shows  part  of  the  Mach  number  sort  of  a  linear  theory  summary. 
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Linear  Theory  Coefficient  Plot 


TABLE  10.  LINEAR  THEORY  SUMMARY  OUTPUT 


*•«  « 

*4  a 


a 

o  o 

U  to 


e 

a  • 

u  u 


to  -» 

V*  fl 

a  u 
•  •— 
•4 


o  o  o  o  o  o 


r*  r*  r*  r»  r»  r» 


o 

m  • 


a 

o  — 
-to  •  r» 

•  to*  —  ■ 

6  At  N 

•  to  a  •* 

O  45  >  • 

•  V 

to  — 


6  W 

O  •  I 

•to  \ 


to  •  O  • 

u  ji  "x  n 

3U  — 


a  >i 

o  o  — 

a  to 

to  i-N 


U 

•  to  —  — 
to  Ih 


0  >i 

o  o  — 
•toe  ■ 
to 


to  ton  • 

a 

to  to  —  — 
to 


I  I 


e 


1  « 

ai  i 

o 

o 

o 

o 

a 

a 

i 

i  o 

1 

0  1 

o 

o 

o 

o 

a 

o 

i 

1  H 

a 

1  to  >• 

U  1 

o 

o 

o 

o 

a 

o 

i 

l  to 

e 

— 

a  a  a 

a 

1  to  m 

1  1 

1  to 

•  • 

• 

• 

• 

1  m 

•to 

■ 

a  r*  a  o 

a 

1  • 

fl  1 

o 

o 

o 

a 

a 

a 

1  V  to 

1  to 

0. 

— 

N 

a  a  a  a 

a 

to 

1  R 

V  1 

a 

o 

o 

o 

a 

a 

1  a  to 

a  a 

a 

a 

a 

1  u 

■ 

a 

a 

a  a  a  a 

a 

? 

1  H 

**  1 

• 

• 

• 

• 

• 

• 

1  O  A 

a  a 

a 

a 

tv 

1  • 

to  *4 

a 

o  o  o  o 

a 

a. 

1 

e 

o 

e 

o 

a 

a 

1  ■  ■ 

At  a 

a 

a 

a 

1  to 

O 

•  •  •  • 

• 

to 

O 

to 

1  « 

1  >%9 

a  a 

o 

a 

a 

l  to* 

a  a  a  a 

a 

a 

• 

1  -to 

«—•  1 

1  •  M 

a  a 

a 

O 

a 

1 

1  1  1 

o 

to 

A 

1  to 

N  I 

•* 

•* 

a 

a 

a 

a 

i  m 

a 

a 

a 

1 

• 

m 

1  to 

fl  1 

• 

• 

• 

• 

• 

• 

i 

1 

>*a 

a 

l  •  M 

U  1 

i 

1 

to 

fl 

■ 

i  a 

1  1 

•* 

•* 

a 

a 

a 

a 

i 

1 

totototo 

a 

to 

a 

1  H 

fl  1 

•to 

•to 

a 

At 

a 

i 

1  fl 

a 

a  a  a  a 

At 

«  to 

JZ 

1 

o*  l 

At 

•to 

•to 

a 

a 

a 

1  «W 

l  o 

c 

— 

a  r-  a  a 

o 

• 

u 

1 

| 

1  o  — 

1  -H 

•to 

fl 

a 

a 

to 

1 

1  *0  u 

At  At 

a 

a 

a 

1  to 

ftto  X 

a  a  a  a 

a 

*» 

a 

ae 

1 

1  Tl  C  • 

a  o 

a 

a 

o 

1  « 

fl  J 

a 

•  •  •  • 

• 

o 

1  to 

1 

1  •  9  to 

•  • 

• 

• 

• 

1  to 

to 

«-• 

• 

>*  to 

>* 

1  **t 

M  1 

a 

a 

a 

a 

a 

a 

1  •  o  \» 

a  a 

a 

a 

a 

1  ? 

o 

1  1  1 

to 

a 

A 

1  to 

■  1 

• 

• 

• 

• 

• 

• 

i  a«*  ■ 

a  a 

a 

a 

a 

i  m 

a 

O 

1  to 

V  | 

o 

o 

o 

o 

a 

a 

i  a  — 

a  a 

a 

a 

a 

i 

• 

to 

to 

1  •  9* 

1  1 

#* 

a 

a 

a 

a 

a 

i 

i 

JS 

e 

to 

1  B 

fl  1 

At 

•to 

•to 

a 

a 

a 

i 

i 

H 

e 

o 

1  M 

e*  i 

At 

•to 

•to 

a 

a 

a 

i 

1  c 

«n 

1 

—•  i 

l 

i  to 

to 

1 

i  >%  — 

o  o 

o 

o 

o 

1  -H 

*• 

« 

1 

1  to  a 

a  a 

a 

a 

a 

1  « 

o 

— 

•«• 

•  •Art 

a 

• 

1 

1  to-H  fl 

a  a 

a 

a 

a 

1  to 

to 

a 

a 

a  rt  «a  a 

At 

e 

l  to 

—  i 

O 

o 

o 

o 

o 

a 

1  a  •  u 

a  a 

a 

a 

a 

1  • 

u  a 

• 

•  •  •  • 

• 

•H 

1  H-H 

N  1 

r» 

r* 

r* 

r* 

r* 

r* 

I4BS 

a  a 

a 

a 

a 

1  U 

to 

to 

a 

ai  a  *1  a 

At 

J 

1  to  to 

fl  1 

•to 

•to 

•to 

a 

a 

a 

1  •  fl 

a  a 

a 

a 

a 

1  •  > 

•a 

a 

I  -to  to 

O  1 

• 

• 

• 

• 

• 

• 

i  a  a 

a  a 

o 

o 

o 

1  to 

1  M  • 

1  1 

m 

to 

a 

a 

a 

a 

i  «• 

•  • 

• 

• 

• 

i  a 

i  <  e 

1  H 

fl  1 
to  1 

to* 

tol 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

1 

1 

a  a 

a 

a 

a 

1 

l 

1 

1 

• 

1 

1 

i 

i  e 

a  a  a  a 

a 

1 

1  Tl  to 

i  • 

to 

v*  a  a  v* 

At 

1 

*—  1 

a 

a 

a 

a 

m 

a 

1  to  U 

1  to 

O 

— 

•  •  •  • 

• 

1 

«  1 

r* 

r* 

r- 

a 

a 

a 

i  >  a 

At  a 

a 

a 

a 

1  to 

to 

a 

a 

a  a  a  a 

a 

1  « 

fl  1 

• 

• 

• 

• 

• 

• 

1  to  e 

•  • 

• 

• 

• 

l  to 

u  a 

• 

1  « 

to  l 

a 

a 

a 

a 

a 

a 

•  •  to  — 

a  a 

a 

a 

a 

l  to 

•  a  ti 

1  to 

to  | 

a 

to* 

a 

a 

a 

a 

1  «  «  ■ 

a  a 

a 

a 

a 

1  9 

> 

a 

•  r. 

V  1 

rt 

to 

a 

a 

a 

a 

1  J»to  - 

a  a 

a 

a 

i  r. 

i 

w  I 

i  o  a 

i 

i 

i 

I 

o 

a 

Ifl  VMfl 

« 

l  • 

i 

i 

a 

• 

•  •  •  • 

• 

1  H  to 

l  to 

i 

a 

• 

*i  a  At  a 

At 

i  to  • 

ai 

a 

a 

a 

a 

a 

i  «m  a 

I 

a 

a 

a  •  h 

At 

1  to  to 

M  1 

*• 

a 

a 

a 

a 

a 

l  •  o 

• 

• 

1  u  • 

■  1 

tol 

a 

a 

a 

a 

a 

I  a 

i 

Tl 

1  •  1 

U  1 

to* 

a 

a 

a 

a 

a 

1  •  to 

a  At 

a 

a 

a 

I 

to 

a 

a  a  a  a 

At 

1  n  to 

«—  1 

• 

• 

• 

• 

• 

• 

1  o  to 

a  a 

a 

a 

a 

1  A 

• 

•  Hrt  A 

r* 

•  O  to 

fi 

« 

a 

a 

a 

a 

1  *  to 

1  u 

a 

•f  rt  rt  rt 

a 

i  to  a 

i  a 

1  to 

■ 

.  •  •  • 

• 

i  a 
i 

I 

I 

• 

i 

9 

a  a  a  a  o 

i 

i 

• 

a 

r» 

a 

a 

a 

I 

i 

a  a 

r* 

•A 

a 

i 

i 

a  a  a  a 

rt 

i 

r* 

a 

r* 

a 

a 

a 

I 

r»  a 

r» 

a 

a 

i 

AaAaa 

1  to 

to 

a 

a 

a 

a 

a 

1  to 

a  a 

a 

a 

a 

i 

to 

At  a  a  At 

m 

1  to  • 

o 

o 

a 

a 

a 

a 

1  to  • 

a  a 

a 

a 

a 

l  to 

to 

a  a  a  • 

a 

1  •  A 

r* 

a 

r» 

r* 

r* 

I  •  A 

a  a 

a 

a 

a 

l  • 

A 

A  A  A  A 

r* 

1  A  fl 

«• 

a 

a 

a 

a 

M 

1  A  fl 

a  a 

a 

0 

1  A 

fl 

a  a  m  m 

a 

I  *4  a 

i  a  a 

i  a 

9 

a 

i 

a 

m 

a 

m 

a 

m 

a 

m 

a 

m 

a 

a 

a 

i 

a  a 

m  m 

a 

a 

a 

a 

a 

m 

i 

i 

K 

a  a  a  a 

a  a  a  a 

a 

a 

« 


83 


TABLE  10.  LINEAR  THEORY  SUMMARY  OUTPUT  (CONCLUDED) 
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SECTION  VII 


SIX  DEGREE  OF  FREEDOM  (6  DOF) 


1.  SIX  DOF  INPUT 

The  6  DOF  reduction  is  split  into  separate  input,  analysis,  and  output 
programs.  The  input  program  creates  a  named  data  file  for  execution  by  the 
analysis.  The  analysis  will  always  run  the  latest  data  configuration  from  the 
input  program. 

The  first  step  in  setting  up  the  input  file  for  the  6  DOF  analysis  program 
is  to  select  the  number  of  shots  to  be  reduced  together.  Up  to  five  shots  may 
be  reduced  in  a  multiple  fit: 

Eater  number  of  shots  in  6  DOF  analysis: 

Next  the  desired  shots  are  identified: 

Key  in  10  digit  6  DOF  projectile  title  (or  LIST) 

:  BS87040818 

The  LIST  option  will  display  all  of  the  shot  numbers  available  for  a 
specified  shot  group.  A  negative  shot  group  number  will  display  all  shots. 

An  automatic  fit  option  is  available  on  request.  This  option  will  allow 
logic  within  the  program  to  turn  on  fit  flags  based  on  the  reduction  and 
projectile  configuration: 

Do  you  want  to  exercise  the  automatic  fit  option  (Y/N)?  N 

If  the  selected  shot(s)  have  been  reduced  before,  the  initial  estimates 
for  this  reduction  may  be  those  used  for  the  previous  reduction.  If  this 
option  is  selected,  the  accuracy  of  the  previous  reduction  will  be  displayed 
and  the  user  prompted  on  whether  or  not  they  shall  be  used: 

Do  you  want  to  update  initial  estimates  based  on 

previous  6  DOF  reduction  on  specified  shot(s)  (Y/N)?  Y 

Initial  estimate  update  based  on  6  DOF  -  reduction  of: 

Date  8-5EP-87  Time  19:35:18 


Probable  errors  of  fit  to  the  data  were: 


Z 


Y-Z  Theta-Psi  Phi 


(m)  (m) 

0.003751  0.002992 


(deg) 

0.748300 


(deg) 

0. 


♦  i  •  win 


Do  you  want  to  update  (Y/N)?  N 


Three  forms  of  parameters  must  be  defined  for  the  analysis:  the  data 
grouping  and  update  options,  initial  estimates,  and  fit  flags.  Initial 
estimates  and  fit  flags  must  be  defined  for  the  first  fit.  Subsequent  fits 
which  turn  on  additional  coefficients  may  also  be  defined. 

The  grouping  and  update  menu  is: 

Fixed  Plane  Analysis  6 -DOF  Input 

Sectional  Fit  Options  (Single  Fit  Only) 


1  “  Data  points  in  first  section  (HA) .  15 

2  -  Data  points  added  to  each  section  (NB). ..........  10 


6D  Summary  /  Update  Options 


3  -  Group  number  code........ . . .  0 

4  -  6  DOF  -  Summary(l,  2,  or  3) .  3 

5-6  DOF  -  Initial  estimate  update(l,  2,  or  3) .  3 


(1  :  no  update)  (2  :  fit  option  selection) 

(3  :  program  criteria) 

Miscellaneous 

6  -  Integration  time  step  (sec)  .00043000 

Fey  in  item  number  ","  new  value 

Key  in  "LIST  to  relist  "DONE"  to  continue 

:  DONE 

Options  One  and  Two  define  how  many  data  points  within  the  range  will  be 
reduced  together  during  the  fit.  Option  Three  requires  the  user  to  specify 
the  group  code  numbers  for  cataloging  the  6  DOF  results.  Options  Four  and 
Five  signal  the  program  for  updating  the  summary  file.  Criteria  One  causes  no 
update  to  be  made  to  the  summary  file.  Criteria  Two  causes  an  update  based  on 
the  setting  of  two  fit  flags.  Criteria  Three  causes  an  update  to  the  summary 
file  if  a  10  percent  improvement  in  the  probable  error  occurs.  Option  Six 
allows  the  computed  integration  time  step  to  be  changed. 
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The  next  input  option  allows  for  changing  the  initial  estimates  that  were 
passed  to  6  DOF  from  the  linear  theory  reduction  program: 

Do  you  want  to  change  initial  estimates  (Y/N)?  Y_ 

Initial  Estimates  for  Shot  Humber:  BS87040818 

1  -  Cma  :  1.13 

2  -  Cma3  :  O.OOOE+OO 

3  -  Cma5  :  O.OOOE+OO 


• 

32  -  Cnga3  :  O.OOOOE+OO 

33  -  Cmga  :  0.0000E+00 

34  -  Cnga  :  O.OOOOE+OO 

Key  in  item  number  new  value 

Key  in  "LIST  to  relist  "DONE"  to  continue 

:  PONE 

The  above  list  contains  a  partial  selection  of  the  input  parameters.  In 
fact,  the  coefficient  input  is  split  info  two  menus.  An  example  taken  from 
the  second  menu  is: 

1  -  Pitch  Angle  (+  down)  [THETA]  (deg):  -1.87 

2  -  Pitch  Rate  (deg/sec):  -165. 

3  -  Yaw  (+  left)  [PSI]  (deg):  0.930 

4  -  Yaw  Rate  (deg/sec):  -823. 

5  -  Horizontal  Range  [X]  (feet):  11.9 


O.OOOE+OO 
O.OOOE+OO 
O.OOOE+OO 
O.OOOE+OO 
O.OOOE+OO 

Key  in  item  number  ntN  new  value 

Key  in  "LIST  to  relist  "DONE"  to  continue 

:  DONE 


UKU 

CmdA 

CmdB 

CXO— unique 
Cln— uniaue 


r 


k 


Once  the  initial  estimates  are  set  the  initial  fit  options  may  be 
modified: 

Do  you  want  to  change  fit  options  (I/H)?  Y 

Fit  Flags  for  Shot  Number:  BS87040818 
40  MM  HEDP  TUBULAR 

1  -  Caa  :  1 

2  -  Cma3  :  0 

3  -  Cma5  :  0 


• 

32  -  Cnga3  :  0 

33  -  Cmga  :  0 

34  -  Cnga  :  0 

Key  in  item  number  new  value 

Key  in  "LIST  to  relist  "DONE"  to  continue 

:  DONE 

Here  too,  the  initial  fit  flag  menu  takes  two  pages: 

Fit  Flags  for  Shot  Number:  BS87040818 
40  MM  HEDP  TUBULAR 

1  -  Pitch  Angle  (+  down)  [THETA]  :  1 

2  -  Pitch  Rate  [THETA-DOT]  :  1 

3  -  Yaw  (+  left)  [PSI]  :  1 

4  -  Yaw  Rate  [PSI-DOT]  :  1 

5  -  Horizontal  Range  [X]  :  1 


• 

17  -  CYgO  :  0 

18  -  CngO  :  0 

19  -  CngO  :  0 

20  -  CXO— unique  :  0 

21  -  Cip-unique  :  0 


Key  in  item  number  n,n  new  value 

Key  in  "LIST  to  relist  "DONE"  to  continue 

:  DONE 

Once  the  initial  fit  flags  are  set  the  user  may  define  up  to  10  subsequent 
sets  of  fit  flags.  The  menu  for  subsequent  fits  is  different  from  the  menu 
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for  the  initial  fit.  Appendix  B  contains  a  complete  list  of  the  coefficients 
and  the  associated  line  numbers.  Here  all  options  appear  on  one  page  and  not 
all  options  are  available: 

*  "  Do  you  want  subsequent  fit  options  (Y/N)? 


Fit  Option  #  2 

1  -  Caa  :  1 

2  -  Cma3  :  0 

3  -  Cma5  :  0 

4  -  Cnq  :  1 

5  -  Caq2  :  0 


32  —  Cnga3  :  0 

33  -  Cnga  :  0 

34  —  Cnga  :  0 

35  —  Caa- unique  :  0 

36  —  Inertial  For*  Factor  [IX/IY]  :  1 

37  —  Cl-DEL  fin  cant  roll  oooent  :  0 

38  -  CHdA  :  0 

39  -  CNdB  :  0 

40  -  CadA  :  0 

41  -  CadB  :  0 

42  —  CXO-unique  :  0 

43  -  Clp-unique  :  0 


Key  in  item  number  new  value 

Key  in  "LIST  to  relist  "DONE"  to  continue 

:  DONE 

Do  you  want  subsequent  fit  options  (Y/N)?  N 

The  input  program  is  terminated  as  soon  as  no  further  subsequent  fits  are 
requested. 

2.  SIX  DOF  ANALYSIS 

The  6  DOF  Analysis  program  may  be  run  in  either  interactive  or  batch  mode. 
When  running  in  a  batch  mode,  the  user  will  be  prompted  to  enter  a  unique  file 
name.  The  input  data  file  will  be  copied  to  a  file  with  this  name.  This  is 
necessary  to  maintain  unique  file  names  when  multiple  processes  are  running. 
When  running  in  the  interactive  mode,  the  program  will  run  the  data  as  last 
configured  by  the  input  program.  The  user  will  see  the  input  data  echoed, 

^  parameter  values  input  during  the  analysis,  and  an  analysis  summary.  The 

program  may  be  aborted  out  of  the  interactive  mode  by  typing  <CTRL  C>.  The 
user  may  then  restart  the  analysis  in  a  batch  mode. 


* 
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The  interactive  analysis  mode  gives  the  user  the  capability  of  modifying 
the  fit  during  the  analysis.  Using  the  initial  estimates  and  fit  options 
selected  in  the  input  program,  a  reduction  is  performed  and  the  probable 
errors  printed  out: 

Probable  Errors  of  Fit  —  Fit  Option:  1 


Pitch-Taw 

Roll 

Travel  [X] 

Swerve  [Y-Z] 

Fit 

(deg) 

(deg) 

(feet) 

(feet) 

1 

0.40402 

0.00000 

0.00531 

0.00893 

2 

0.33277 

0.00000 

0.01081 

0.00549 

3 

0.32895 

0.00000 

0.00162 

0.00534 

4 

0.32893 

0.00000 

0.00161 

0.00534 

Following  the 

probable  errors , 

the  values  of 

the  coefficients  used  for 

each  fit  will  be  printed: 

S (unary  of  Aerodynamic  Coefficients  During  Fit 


1 

Cma 

1.26880 

1.27731 

1.27708 

1.27707 

4 

Cmq 

-0.11800 

0.28712 

0.27919 

0.27675 

6 

Cnpa 

-0.23880 

-0.42858 

-0.43711 

-0.43760 

17 

CXO 

0.16070 

0.15089 

0.15072 

0.15088 

20 

CNa 

1.90340 

1.99709 

2.00090 

2.00133 

32 

Clp 

-0.02090 

-0.01120 

-0.01186 

-0.01186 

The  coefficients  should  be  examined  to  see  if  a  solution  is  being 
converged  on.  If  the  values  from  fit  to  fit  are  oscillating,  the  parameter 
should  be  fixed. 

After  this  table  is  printed,  another  table  of  the  fit  parameters  and  their 
associated  probable  errors  will  be  displayed: 

Final  Coefficients  —  Fit  Option:  1 
At  average  Mach  number  unless  noted 


Parameter 

Probable  Error 

Fit  Flag 

1  Cma 

1.27707 

0.01158 

1  :  201 

4  Cmq 

0.27675 

0.29260 

1  :  204 

6  Cnpa 

-0.43760 

0.09651 

1  :  206 

7  Cnpa3 

5.54210 

0.00000 

0  :  207 

17  CXO 

0.15088 

0.00208 

1  :  217 

18  CXa2 

2.30000 

0.00000 

0  :  218 

20  CNa 

2.00133 

0.03661 

1  :  220 
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Final  Coefficients  —  Fit  Option:  1  -  (Continued) 

23  CTpa  :  >1.00000  0.00000  0  :  223 

32  Clp  :  -0.01186  0.00525  1  :  232 

146  Fit  section  :  10 

147  Section  increment  :  5 

Total  data  points  :  26 

Key  in  itea  number  ","  new  value 
"LIST"  to  relist  "ABORT"  to  abort 

"DONE"  to  increment  fit  "HELPLESS"  to  end  HELP 
"RERUN"  for  same  fit  length 


At  this  point  in  the  program,  parameters  may  be  changed  again.  Fit 
options  may  also  be  turned  on  or  off.  The  numbers  along  the  left  side  allow 
the  user  to  modify  a  parameter.  The  numbers  along  the  right  allow  the  user  to 
change  fit  flags.  All  non-zero  parameters  will  be  displayed.  Values  or  fit 
flags  other  than  those  displayed  may  be  changed  but  the  user  must  either  know 
the  parameter  number  or  refer  to  the  list  in  Appendix  B.  Note  the  numbers  in 
the  right  column  are  the  same  as  the  left  column  with  200  added. 

The  RERUN  option  will  allow  the  analysis  to  be  repeated  for  the  same 
number  of  tunnel  stations  with  the  modified  parameters  and  rlags.  The  DONE 
option  will  also  use  the  modified  parameters  and  flags  but  will  increment  the 
number  of  stations  by  the  summing  section  increment.  The  ABORT  option  will 
terminate  the  analysis  program.  The  HELPLESS  option  will  allow  the  analysis 
to  continue  from  this  point  without  displaying  this  change  menu  again.  The 
program  will  complete  this  fit  and  any  subsequent  fits  as  defined  by  the  input 
program. 

The  change  menu  will  only  appear  during  the  first  fit  of  the  data.  Fit 
options  during  subsequent  fits  come  from  the  input  data  file. 

3.  SIX  DOF  OUTPUT 

The  6  DOF  output  program  allows  for  the  creation  of  summary  tables,  motion 
plots,  and  coefficient  plots.  The  output  menu  is: 

AERODAS  6  DOF  Output  Segment 


1  Summary  Tables  (Terminal) 

2  Summary  Tables  (Llneprinter) 

3  Motion  Plots 

4  Coefficient  Plots 

5  Exit  Output  Program 


If  summary  tables  are  requested,  the  shot  group  to  be  output  must  be 
specified.  The  summary  file  will  be  sorted  for  the  single  and  multiple 
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reductions  for  this  group  number.  Ihv  user  may  select,  from  the  shots,  which 
ones  to  Include  in  the  summary,  liable  ii  shows  a  summary  of  40mm,  spin 
stabilized,  shot  reductions. 

Coefficient  plots  are  also  created  by  shot  group  numbers.  For  the  shots 
selected  from  the  desired  shot  group,  the  following  menu  of  plots  is 
available: 

Summary  Segment  -  Six  Degree  of  Freedom  Reductions 

SHOT  GROUP  1  —  Sample  Output 

vs  Mach  Number  vs  Angle  of  Attack 


1  -  cxo 

CX  -  11 

2  -  CNa 

CN  -  12 

3  —  Cma 

Cm  -  13 

4  -  Cnpa 

Cap  “  14 

5  -  Cmqa 

Cmq  -  15 

6  -  Clp 

7  -  CXa2 

8  —  Cma  3 

9  —  Cnpa3 

99  -  Done 

Plots  versus  Mach  number  are  for  the  single  shot  reductions.  Plots  versus 
angle  of  attack  are  for  multiple  shot  reductions.  Figures  20  and  21  show  drag 
coefficients  in  both  plot  modes. 

When  a  motion  plot  is  to  be  created,  the  program  reads  the  coefficients 
for  the  specified  single  or  multiple  shot  reduction.  These  coefficients  are 
then  displayed  for  modification  by  the  user.  Each  of  the  shots  will  be 
integrated  through  its  trajectory.  Probable  errors  for  each  shot  in  the 
reduction  will  be  displayed  prior  to  the  display  of  the  following  plot  menu: 

Plot  Segment  -  Six  Degree  of  Freedom  Reductions 

SHOT  GROUP  1  —  Sample  Output 

Plot  #  Description  Plot  #  Description 

1  Y  vs  X  8  Residual  in  X  vs  X 

2  Z  vs  X  9  Residual  in  Y  vs  X 

3  Phi  vs  X  10  Residual  in  Z  vs  X 

4  Theta  vs  X  11  Residual  in  Phi  vs  X 

5  Psi  vs  X  12  Residual  in  Theta  vs  X 

6  Alpha  vs  X  13  Residual  in  Psi  vs  X 

7  Theta  vs  Psi  14  Velocity  vs  X 

19  Theta  &  Psi  vs  X 

20  Change  Shot  Number 

21  Print  tabulated  data  on  exit 

99  Done 
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TABLE  11.  6  DOF  SUMMARY 
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TABLE  11.  6  DOF  SUMMARY  (CONCLUDED) 
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DOF  Reduction  -  Coefficient  Output 


C_>  X 
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Figure  20.  6  DOF  Coefficient  Versus  Angle  of  Attack 


DOF  Reduction  -  Coefficient  Output 


OXO 
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Figure  21.  6  DOF  Coefficient  Versus  Mach  Number 


From  the  above  menu  either  motion  or  residuals  may  be  plotted  against 
the  travel  down  the  range.  Figure  22  shows  a  motion  plot.  Figure  23  shows  a 
residuals  plot.  For  a  multiple  shot  reduction,  the  data  from  each  shot  is 
»  available.  The  program  will  default  to  the  first  shot  in  the  reduction  for 

plotting.  Option  Twenty  enables  the  user  to  select  other  shots  for  plotting. 
Option  Twenty  One  enables  tabulated  data  from  the  motion  integration  to  be 
printed  upon  exit  from  the  program. 
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BS87032797  29-JUL-87  12:12:24  40  MM  HEDP  TUBLAR 

6  DOF  Reduction  -  Motion  Plot 


Q.  (/)  H 


HXUJI-C 
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Figure  22.  6  DOF  Angular  Motion  Fit 


BS87032797  29-JUL-87  12:12:24  40  MM  HEDP  TUBLAR 

6  DOF  Reduction  -  Motion  Plots 


Figure  23.  6  DOF  Data  Residual 
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APPENDIX  A 


TECHNICAL  BACKGROUND  AND  EQUATIONS  OF  MOTION 


APPENDIX  A 


TECHNICAL  BACKGROUND 
AND 

EQUATIONS  OF  MOTION 


This  section  provides  the  equations  of  motion  contained  in  the  linear  theory 
and  6DOF  programs  in  ARFDAS.  These  mathematical  models  are  directly 
correlated  to  the  spark  range  data  resulting  in  a  determination  of  the 
aerodynamic  coefficients  and  coefficient  derivatives  that  provide  the  best  fit 
for  the  experimental  data. 

Correlation  of  the  data  to  the  linear  theory  equations  of  motion  is  done  using 
a  standard  least  squares  method.  This  analysis  provides  for  preliminary 
analysis  and  screening  of  the  data  prior  to  the  6  DOF.  Initial  estimates  for 
the  initial  conditions  and  coefficients  are  obtained  based  on  this  analysis. 

A  more  detailed  description  and  derivation,  of  the  equations,  is  contained  in 
Reference  A-1. 

Correlation  of  the  data  to  the  6  DOF  equations  of  motion  is  performed  using 
the  Maximum  Likelihood  Method.  There  are  two  forms  of  the  dynamic  equations 
available  in  ARFDAS.  Selection  of  the  set  of  equations  depends  on  the 
configuration  being  analyzed.  For  a  configuration  with  both  physical  and 
aerodynamic  symmetry,  the  dynamic  equations  derived  in  the  fixed  plane 
(nonrotating)  coordinate  system  are  used  (References  A-2  and  A-3).  When 
either  the  physical  or  aerodynamic  model  contain  asymmetries,  the  dynamic 
equations  derived  in  the  body  fixed  (rotating)  coordinate  system  are  used 
(Reference  A-i|).  It  should  be  noted  that  both  6  DOF  programs  contain  a 
generalized  aerodynamic  model.  The  significant  coefficients  and  derivatives 
to  be  determined  depend  on  the  type  of  configuration  being  analyzed. 

Solution  of  the  dynamic  equations  is  performed  using  a  fourth  order  Runge- 
Kutta  integration.  Useful  transformation  identities  between  various 
coordinate  systems  have  been  included  in  this  appendix.  The  sign  convention 
for  the  fixed  plane  coordinate  system  is  illustrated  in  Figures  A-1  and  A-2. 
Typical  forms  of  the  induced  aerodynamic  coefficients  for  a  four-fin  missile 
are  illustrated  in  Figures  A-3  and  A-4. 
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DOWN  RANGE 


LOOKING  DOWN  RANGE 


SIDE  VIEW 

POSITIVE  DIRECTIONS  SHOWN 


Figure  A-l.  Coordinate  System  Fixed  Plane  (Positions) 
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LOOKING  DOWN 
RANGE 


Figure  A-2.  Coordinate  System  Fixed  Plane  (Rotation) 
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Figure  A-3.  Roll  Angle  Orientations 
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Figure  A-4.  Illustrations  of  Induced  Effects  (4  Fin  Missile) 
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1.  NOMENCLATURE 


Physical  Properties 


A 

Area  (ird2/4) 

ft2 

d 

Diameter 

ft 

I 

X 

Axial  moment  of  inertia 

slug-ft2 

*» 

Transverse  moment  of  inertia 
(y-axis) 

2 

slug-ft 

Transverse  moment  of  inertia 
(z-axis) 

slug-ft2 

m 

Mass 

slugs 

CG 

Center  of  gravity 

1  from  nose 

L 

Projectile  length 

ft 

m 

c 

Mass  of  ball  rotor 

slugs 

rx 

Distance  from  CG  to  CG  of  ball  rotor 

ft 

ryz 

Radial  clearance  of  ball  rotor 

ft 

Subscripts 

r 

Reference  value 

Range  Properties 

• 

P 

Atmospheric  pressure 

millibars 

T 

Temperature 

Degrees  celdus 

RH 

Relative  humidity 

Z/100 

R 

n 

Reynolds  number 

— 

M 

Mach  number 

— 

a 

Speed  of  sound 

ft/sec 

p 

Air  density 

3 

slugs/ft 

g 

Gravity 

32.134  ft/sec2 
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nT 


Range  Properties  (continued) 


30.5  degrees 
126.0*  degrees 

e 

Subscripts 

r  Reference  value 


Range  latitude 

6  Range  aziouth 

R 

w  Earth  rotation  rate 


Linear  Theory 

dl,2,3,4 
rl»2, 3,4 


Constants  for  drag  equation 
Constants  for  roll  equation 
Amplitude  nutation  vector,  degrees 
Amplitude  precession  vector,  degrees 
Amplitude  trim  vector,  degrees 
Amplitude  yaw  of  repose,  degrees 
Orientation  nutation  vector,  degrees 
Orientation  precession  vector,  degrees 
Orientation  trim  vector,  degrees 
Angular  velocity  nutation  vector,  degrees/ft 
Angular  velocity  precession  vector,  degrees/ft 
Angular  velocity  trim  vector,  degrees/ft 

2 

Angular  acceleration  nutation  vector,  degrees/ ft 

2 

Angular  acceleration  precession  vector,  degrees/ft 
Damping  factor  nutation  vector,  1/ft 
Damping  factor  precession  vector,  1/ft 
Amplitude  nutation  swerve  vector,  ft 
Amplitude  precession  swerve  vector,  ft 
Orientation  nutation  swerve  vector,  degrees 
Orientation  precession  swerve  vector,  degrees 
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Linear  Theory  (continued) 


* 


9 


zs 


1,2,3 
ySl,2.3 


Subscripts 


m 


Constants  for  swerve  equation  (2 -plane) 

Constants  for  swerve  equation  (y-plane) 

Distance  from  reference  station,  ft 

Distance  from  reference  station  to  midpoint  of  fit,  ft 
Midpoint  of  fit,  ft 

Evaluated. at  midpoint  of  fit 


* 


6 DOF  -  Fixed  Plane 
x.y,z 

u,v,w 

P.q,r 


a  , 

a  , 

a 

ex* 

cy* 

cz 

a  , 

a  , 

a 

cu* 

cv* 

cw 

a 

e 


Y 

V 


F 

z 


L,  M,  N 


Missile  coordinates  (fixed  plane)  ft 
Missile  velocities  (non-rolling)  ft/sec 
Missile  angular  velocities  (non-rolling)  rad/sec 
Missile  orientation  (fixed  plane)  radians 

2 

Components  of  coriolis  acceleration  (fixed  plane),  ft/sec 

2 

Components  of  coriolis  acceleration  (non-rolling),  ft/sec 
Total  angle  of  attack,  radians 
Sine  of  the  total  angle  of  attack 
Aerodynamic  roll  angle,  radians 
Total  missile  velocity,  ft/sec 
Force  components  (non-rolling) ,  11 
Moment  components  (non-rolling),  1 


Superscript 


First  derivative  with  respect  to  time 
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£52£  --Body- Fixed 


X,  y,  2 

V  V  wb 

V  qb*  rb 

9,*,4 

acx»  acy*  acz 
acub*  acvb*  acwb 


Missile  coordinates  (fixed  plane),  ft 
Missile  velocities  (rolling),  ft/sec 
Missile  angular  velocities  (rolling),  rad/sec 
Missile  angular  orientation  (fixed  plane),  radians 

2 

Components  of  corlolis  acceleration  (fixed  plane),  ft/sec 

2 

Components  of  corlolis  acceleration  (rolling),  ft/sec 


a 


c 

Y 

Fxb*  Fyb*  Fzb 

V  V  *b 


Total  angle  of  attack,  radians 
Sine  of  the  total  angle  of  attack 
Aerodynamic  roll  angle,  radians 
Force  components  (rolling),  lbs 
Moment  components  (rolling),  ft-lbs 


Superscripts 


First  derivative  with  respect  to  time 


Aerodynamic  Coefficients 


s 

Sc 

CN 

s 

Sy 

CNy 

Sy 

S 


C 


Z 


So 


Drag  coefficient 
Axial  force  coefficient 
Normal  force  coefficient 
Magnus  force  coefficient 
Induced  axial  force  coefficient 
Induced  normal  force  coefficient 
Induced  yav  force  coefficient 
Yaw  force  coefficients 
Normal  force  coefficient 
Trim  yaw  force  coefficient 
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Aerodynamic  Coefficients  (continued) 


♦ 


9 


Zo 

SmV  Stt^B 

s 

C16 

CiP 

CiY 

Cm 


mq 


nr 


my 


nY 


mo 


no 


Co66A*  Cm66B 


Let: 


Trim  normal  force  coefficient 

Trim  force  coefficient  components  (fixed  plane) 

Magnus  force  coefficient 

Roll  moment  coefficient 

Roll  deceleration  coefficient 

Induced  roll  coefficient 

Pitching  moment  coefficient 

Yaw  moment  coefficient 

Pitching  damping  coefficient 

Yaw  damping  coefficient 

Induced  pitching  moment  coefficient 

Induced  yaw  moment  coefficient 

Trim  pitch  moment  coefficient 

Trim  yaw  moment  coefficient 

Trim  moment  coefficient  components  (fixed  plane) 


a 

0 

«,  e 


VT 


Vv 


■  sin 


-1  (y2  t  w2) 


Subscripts 

o^  derivative  with  respect  to  a* 
8^  derivative  with  respect  to  8* 
derivative  with  respect  to  e* 
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2.  LINEAR  THEORY 


Linear  Theory  has  the  following  features: 

o  Self  start  (requires  only  the  precession  frequency  to  be  reasonably 
guessed) 

o  Elimination  of  erroneous  data  points 


The  equations  shown  below  define  the  notion  for  drag,  roll,  yaw  and 
swerve. 


Drag  Analysis 


'  ■  41  +  d2  xr  +  i3  xr  +  d4  *r 


Roll  Analysis 


2  3 

r,  +  r«x  +  r,x  +  r.x 
1  2  r  3  r  4  r 


Yaw  Analysis  2 

X  x  x 

®m  ’  «  1  *“*  <*1  +  *1  xr  +  *1  ~T  > 

.  2 

X2xr  .  -  *  r 

+  ^  e  cos  (#2  +  42  xr  +  *2  ~  * 

+  K3  cos  (/J  f3) 

X,x  .  ..  x2 

*»  -*!•  •!•  (*!  +  *i  *r  +  ~T  > 

X2xr  .  ~  x2r 

+  K2  e  *  rsin  (42  +  42  xr  +  *2  “F  > 

+  K3  sin  (/*  43) 

+  K4 
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..  x 


y  *  y8!  +  ys2xr  +  sie  rsin  (*si  +  *1  xr  +  *1  "T* > 


.  «> 

AX  x 

+  S2e  2  rsin  (*s2  +  *2  \  +  *2  T~  ) 


2  XlXr  •  **  x2 

z  -  zsx  +  zs2xr  +  zs3  xr  +  S^e  cos  (+sl  +  ^  xr  +  ^  ^_r  ) 


X  x  x 

+  S2e  2  rcos  <^S2  +  J2  xr  +  J2  ) 


The  reference  station  for  the  analysis  is  Station  1 
Let  xf  ■  xt  -  Xj^  (distance  from  reference) 


Definition  of  Terms 

The  unknown  coefficients  which  are  solved  are: 
d1#  d2,  d3,  d^  -  drag  analysis 

rl*  r2*  r3*  rA  “  ro11  analysia 

•  •  **  '* 

^2*  ^1*  ^2*  ^2 9  ^l9  ^2 9  9  ^2 9 


K3.  *3 


-  yaw  analysis 


sl»  s2»  +si»  *S2  »  ysl’  ys2f  zal'  zs2*  za3  “  8werve 


^3  is  set  equal  to  the  experimental  roll  position 


<»1  +  *1  ».  +  *2  +  *2  V  « 

<*1  +  *i  V<*2  +  +2  xm>  V. 


2  (yew  of  repose) 
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0  ■ 
ra 


6  - 


£ 

X 


m 


y  +  1 


In  order  to  compute  the  aerodynamic  coefficients,  the  determined  constants 
are  shifted  to  the  midrange  distance. 


m 


(xlast  "  xfirst*  ’  2 


x  -  x  (distance  to  be  shifted) 
iq  r 


♦lm  “ 

♦l  +  *1  xs 

*2.  “ 

*2  +  *2  *3 

•• 

♦lm  “ 

•  • 

♦l 

♦2m  * 

•  • 

*2 

*!«- 

Alxs 

V 

K2m“ 

Vs 

K2e  1  8 

S2m  “ 

s  Z2*3 

S2e 

d2mm 

d2  +  2d3  *  xs 

+  3d 

d3m“ 

d3  +  3  •  dA  • 

xg 

r2m  “ 

r2  +  2r3  •  Xg 

+  3r 

r3m 

r3  +  3  •  r4  • 

xs 
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CD  *  K  <d3„'d2»> 


5  .-I 


Ctp  *  2d  <r3«,/r2»  *  d3»/d2»>  ( ST  >  -  ro11  fit 


c,p  *  2d  »*ln  +  *2m)  7  -  d3„/d2„>  (  llr-)’ -  «* 


UJ 


.  .  .  2  "  3m  2ar  V  161 

♦la  +  *2. 


C  -  = 
ma 


8IV  *lm  ;  *2m 
trpd8 


“  (X1  +  X2  + 


pAV  16_i 


-) 


trpd 


Cnpa  *  FX  -  !r  a  +  +  C*a  (1  -  '»  Ix 


.2 

md  t 


(*lm  +  *2m)/(*lm  “  ^m5 


Sr  S  + 


^2m  d  S2m 
^  ^m  5 


No  attempt  is  made  during  the  linear  theory  analysis  to  compute  C£6  or 
as  these  will  be  determined  in  the  6  DOF  analysis.  These  coefficients 
are  not  necessary  as  estimates  for  the  6  DOF  analysis. 
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3.  SIX  DEGREE  OF  FREEDOM 


The  data  analysis  systm  currently  consists  of  two  6  DOF  programs. 


o  Fixed  plane  6  DOF  for  symmetric  missiles  or  projectiles  including 
induced  forces  and  moments. 


o  Body  fixed  6  DOF  for  configurations  such  as  airplanes  or  symmetric 
missiles  with  or  without  mass  asymmetries. 


Both  6  DOF  programs  can  reduce  (analyze)  up  to  five  sets  of  experimental 
data  (shots)  simultaneously  and  determine  the  set  of  aerodynamic 
coefficients  which  best  fit  all  the  data. 


The  equations  of  motion,  force  and  moment  summations,  and  aerodynamic 
coefficients  will  be  discussed  in  this  section  for  each  6  DOF  program 
individually. 
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6  OOF  Fixed  Plane 


The  equations  of  motion  have  been  derived  In  a  fixed  plane  coordinate 
system.  The  X-axis  points  dovnrange,  the  Y-axis  points  to  the  left  looking 
downrange  and  the  Z-axls  points  up.  Coriolis  effects  have  been  accounted  for. 


x  ■  u  cos  8  cos  i -  v  sin  i>  +  w  sin  8  cos 
y  *  u  cos  8  sin  $  +  v  cos  i>  +  w  sin  8  sin 
z  -  -u  sin  8  +  w  cos  8 


a 

cx 


a 


cz 


-2ue  (y  sin  AR  +  z  cos  AR  sin  dR) 

2a>e  (x  sin  AR  -  z  cos  AR  cos  SR) 

2we  (x  cos  AR  sin  6r  +  y  cos  XR  cos  «R) 


cu 


cv 


cw 


*cx  cos  6  coa  *  +  aCy 
-acx,in*  +  .cy 

a  sin  8  cos  *  +  a 
cx  T  cy 


cos  8  sin  ij;  -  a  _  sin  8 
cz 

COS  1>  • 

sin  8  sin  4>  +  a  cos  8 

CZ 


0 

* 

♦ 


q 

r/cos  8 
p  +  rtan  8 


u 

v 

w 


g  sin  8  -  qw+rv-a  +F/m 

CU  X 

-ru  -rwtan  8  -a  +  F  /m 

cv  y' 

-g  cos  8  +  rv  tan  8  +  qu  -  a  +  F  /m 

M  cw  z 
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p  ■  L/Ix 

q  ■  -rSan  0  -  I  /I  n>  +  M/I 

X  7  7 

r  -  +qrtan  0  4*  I  /I  qp  +  N/I 

x  y  y 

Aerodynamic  Model 

The  aerodynamic  forces  and  moments  of  symmetric  spin  stabilized  projectile 
and  fin  stabilized  missiles  are  shown  below  as  they  are  currently  modeled  in  the 
MLMFXPL  program.  For  the  fin  stabilized  case,  the  analysis  Includes  the  Induced 
forces  and  moments.  At  present,  the  coefficients  are  assumed  to  be  functions  of 
Mach  number,  the  sine  of  the  total  angle  of  attack,  and  the  aerodynamic  roll  angle. 

Summation  of  Forces  and  Moments 
Fx  “ 

Fy  q  A  I-CNo  V  +  2 V  Spa  V  +  Sya  V 

•*N«aA  sin  *  '  co*  ♦  1 

« 

Fz  ’  A  t'c8a  V*  27  Spa  V  "  Sra  V 

-Sl5SAco*  ♦  -  Vll  ,ln  ♦  1 

M  ■  <7  Ad  [C  —  +  Si  Q  +  £i  (f  — 

H  lma  V  2V  mq  2V  npa  V 

^nya  7  +  ‘m66AC0S  ♦  “  3m66B  ♦  +  **  > 
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N  -  q  Ad  [-C  7  +  I^C  +  C  £ 

1  no  V  2V  nq  2V  npo  V 


+  W  v  +  c«4{*  *ln  ♦ +  c»i5B  «■  ♦ +  ™Q  i 


Aerodynamic  Coefficient  Expansions 


Sc  -  Sto  -  Ss2  '2  -  Sc4  '4  +  St.  <Mi  -  V  -  ScyS,  ‘2  «a  »  * 


Sfa  ■  Sta  f  Sta3  c  +  S<a5  1  *%azc  “•  S  y  *  Siam  (M1  -  »r> 


Sj66A  "  Sfi6A 


Sj66B  "  Sffi^B 


Sfpa  +  Sfpa, 


"  Cy  ~  c  sin  N  y 


C.  +  C.  -  C  +  c,  (M.  -  M  ) 
Ip  ipa^  1pm  1  r 


cts{  +  ctw2  c  *ln  N  * 


C  —  +  C  —  c*  +  C  —  t*  +  C  (M.  —  M  )  +  C  —  cos  N  v 

TOOL  TOOL ,  mr*  m  '  i  mvn  • 


1  r'  orra. 


+  CL.—  (CC  -  CC  )  +  C  —  c  cos  N  y 
wa  r  my  ae 


c-+  C - c2  +  C  —  c4 

oq  oq  a,  e  mqa,  e 


C  r  +  C  -  c2  +  C _ -  c4  +  C _ -  c6 
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o  4 

C  -  c  -  e  sin  N  y  +  C  -  c  sin  N  y 
nya  nya,  ay  a. 


4  A  "  Co6  4A 


Safi  6B  "  Safi  6B 


The  following  terns  represent  nonplanar  damping  variations. 


RMQ 

RNQ 


•  Coq2K  1  <2^  (^25  ‘  ^  ^  1  +  Cnq2N  [  (27J  t^2)  +  (2V>  iy2)  1 
’  ^mq2N  '  <W>  ‘  #  ^  >  +  Cnr2H  1  #  <£>  +  ^  <$>  1 


The  aerodynamic  roll  angle  is  computed  by  transforming  the  fixed  plane  velo- 
cities  into  the  rolling  body  coordinate  system  as  follows: 


v.  ■  v  cos  4  +  w  sin  ♦ 


w^  ■  -v  sin  4  +  w  cos  t 
Y  -  tan”1  (vb/wb) 


6  DOF  Body  Fixed 

Since  the  aerodynamic  forces  and  moments  are  defined  in  a  body-fixed 
coordinate  system,  the  equations  of  motion  were  derived  with  respect  to  a 
rotating  coordinate  system.  This  coordinate  system  is  defined  with  the 
x-axis  aligned  with  the  longitudinal  axis  of  the  missile  and  points  out 
the  nose,  the  y-axis  points  out  the  left  wing,  and  the  z-axis  points  up 
with  respect  to  the  body.  The  body-fixed  coordinate  system  is  rigidly 
attached  to  the  missile  and  rotates  with  the  missile  about  the  x-axis. 

The  inertial  frame  of  reference  is  the  Earth.  It  is  assumed  that  the 
Earth  is  fixed  in  space  and  flat.  Coriolis  effects  are  included  in  the 
equations. 
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•  N 


X 


cos  0  cos  $  +  (sin  0  sin  4  cos  4  -  cos  4 
+  w.  (sin  0  cos  4  cos  4  +  sin  4 

D 

cos  0  sin  4  +  vb  (sin  0  sin  4  sin  4  +  cos  $ 


sin0 


+  (sin  0  cos  $  sin  4  -  sin  4 
+  vb  cos  0  sin  4 
+  w.  cos  0  cos  4 

D 


CX 


-2ui  (y  sin  X_  +  z  cos  X_  sin  6_) 

C  K  K  K 


cy 


cz 


2u>  (x  sin  X_  -  z  cos  X_  cos  6_) 

C  K  K  K 


2u»e  (x  cos  XR  sin  6R  +  y  cos  XR  cos  6R) 


a 


cub 


a 


cvb 


a  cos  0  cos  4  +  a  cos  0  sin  ip  -  a  sin 
cx  cy  cz 

acx  (sin  0  sin  4  cos  4  -  cos  $  sin  ip) 

4acy  (sin  0  sin  4  sin  ip  +  cos  $  cos  ip) 

+a  (sin  4  cos  0) 
cz 


a 


cwb 


a  (sin  0  cos  0  cos  4  +  sin  4  sin  ip) 

CX 


+acy  (sin  0  cos  4  sin  ip  -  sin  4  cos  ip) 
+acz  (cos  4  cos  0) 


0 

4 

4 


qbcos  4  -  rbsin  4 


(qbsin  4  +  rbcos  4) /cos  0 
Pb  +  tan  0  (qbsin  4  +  rbcos  4) 


rb  Vb  “  %  Wb  “  “cub  +  Fxb/n  +  8  sin  0 

Pb  Wb  "  rb  Ub  "  “cvb  +  Fyb/n  "  8  Sin  ♦  cos  6 

%  “b  ’  pb  vb  '  “cwb  +  Fzb/n  ”  «  cos  *  cos  0 


sin  ip) 
sin  4) 
cos  4) 
cos  4) 
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N.  +  I 
b  a 


V  +  (I« 


V  Vk 


Aerodynamic  Model 

The  aerodynamic  forces  and  momencs  currently  modeled  In  the  body  fixed 
program  are  given  below.  Coefficients  are  assumed  to  be  functions  of  Mach  number, 
the  sine  of  the  pitch  and  yaw  angles,  or  the  total  angle  of  attack  and  the  aero¬ 
dynamic  roll  angle. 


Summation  of  Forces  and  Moments 


-qAc, 


■  q  A  (-  c  -  CVD  — 


-  h  -  h 
SfB  v"  +  °Yya  V 


Pkd  wk 

r  °  .  _Jd.  ) 
i  pa  2V  V  1 


•  i  A  <-  - 


C  Ik  _  c  i 

Za  V  xya  V 


'‘Ypa  2V  V  * 


pbd  - 

" A  d  < ■&-  ctP  +  ct  > 


—  m  _  q,  d  p.  d  v.  v 

•  q  A  d  (C  +  C  +  C  •=—  +  C  _k  +  g  _J>  * 

mo  ma  V  mq  2V  npa  2V  V  +  nya  V  ' 

v.  r.  d  _  p.d  w.  _  w. 

q  A  d  (+  CnQ  -  CRg  ~  +  Cnr  Yv“  +  ^npa  W  V~  +  ^nya  V~ 
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Aerodynamic  Coefficient  Expansions 

2  2 


cx 

’  CXo+  CXa  (V  *  +  CXa2  (V  *  +  Ss2(V  *  +  CXm  (Mi  "  Mr*  +  Sa2  C 

-  So  “ 

'  So 

So  “ 

CZo 

Se 

2 

-Se  -  S03  0  +  Sa  +  CNo3  c2  +  CNY53  COS  N  Y 

So  “ 

w  w  2 

Sa  +  Sa,  0  +  Sa3  0  +  Si  +  Si3  «*  +  S^  «os  N  * 

Sya- 

Sya3  c2  sln  N  Y 

Spa  " 

Spa 

Clp 

Sp 

s  - 

• 

ct« 5  +  c»-,o2  *2  ■ln  * T  +  Cts  <<r> 

So  - 

C 

mo 

C 

no 

C 

no 

C 

ma 

w  w  2 

1  c«. +  c-.2  0  +  c„3  0  +  C^  +  c«;3  ‘2  +  W3  '2  c0‘  M  * 

*  C  . 

raq 

+  S,  (CO  -  CCr) 

„  2 

,b.  2 

■  C  +  C  <TT-)  +  c  -  +  C  --  C 

mq  mqa^  V  m<l  m<la2 
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—  VK  0  2 

Ca  -  -c  a  -  C  a  (“)  +  C  -  +  C  -  €  +  C  -  e  cos  N  Y 

n8  nB  V  ma  ma^  rayo^ 


+  CYg  (CG  -  CGr) 


vb  2 

C  -  C  +.  C  .  fcr^)  +  C  -  +  C  —  e 
nr  nr  nr8~  V  '  raq  mqa_ 


C  «  C  c  sin  N  y 
nyo  nya_ 


C  -  C 
npa  npa 


The  aerodynamic  roll  angle  is  computed  as  follows: 


Y  -  tan"1  (vb/wb) 
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Trana formations  (Fixed  Plane) 


Missile  to  Earth 

x  ■  u  cos  8  cos  ip  -  v  sin  <p  +  w  sin  8  cos  <p 

y  ■  u  cos  8  sin  ip  +  v  cos  xfi  +  w  sin  8  sin  ip 

■  -u  sin  8  +  w  cos  8 

Earth  to  Missile 

u  ■  x  cos  8  cos  ip  +  y  cos  8  sin  ip  -  z  sin  8 

v  *  -x  sin  tp  +  y  cos  ip 

w  ■  x  sin  8  cos  ip  +  y  sin  8  sin  ip  +  z  cos  8 


Transformations  (Body  Fixed) 
Missile  to  Earth 


cos  8  cos  ip 


+  (sin  8  sin  $  cos  ip 

+  (sin  8  cos  $  cos  ip 


+ 


cos  <p  sin  ) 
sin  <p  sin  ip  ) 


y 


cos  8  sin  ip 


+  (sin  8  sin  <p  sin  ip 

+  (sin  8  cos  $  sin  tp 


+ 


COS  <p  cos  ip  ) 
sin  <p  cos  ip  ) 


z 


-u^  sin  8 


+  cos  8  sin  $ 

+  cos  8  cos  $ 
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-  x  cos  6  cos  +  y  cos  8  sin  -  z  sin  6 

*  x  (sin  8  sin  <p  cos  $  -  cos  4  sin  t|>  ) 

+y  (sin  8  sin  4  sin  f  +  cos  4  nos  \(i  ) 

+z  (cos  8  sin  t  ) 


-  x  (sin  8  cos  +  cos  \ji  +  sin  t  sin  ) 
+y  (sin  8  cos  +  sin  \(i  -  sin  <p  cos  ) 

+x  (cos  8  cos  +  ) 
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Transformations 


•) 


t> 


Body  Fixed  to  Fixed  Plane 


u 

"  “b 

V 

"  vb 

coa  4 

wb 

sin  4 

V 

"  vb 

sin  4 

+  wb 

coa  4 

P  -  Pb 

q  -  qb  coa  4  -  rb  sin  4 
r  ■  qb  sin  4  +  *b  coa  4 


Fixed  Plana  to  Body  Fixed 


«b  -  n 

vb  -  v  coa  4  +  w  sin  4 
wb  -  -v  sin  4  +  w  coa  4 


Pb  "  P 

qb  -  q  cos  4  +  r  sin  4 
rb  ■  -q  sin  4  +  r  cos  4 

Direction  Coainaa  to  Fixed  Plane 

0  ■  -  sin  *  (n  ) 

6 

-i  2  2  ^  ^ 

4  -  sin  (ne/(me  +  p#  )  ; 
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Angular  Identities 


sin  a  ■ 

SI 

sin  6  ■ 

n 

sin  0  ■ 

m 

sin  *  - 

to 

tan  0  ■ 

a 


Mlssila  Angles 

/  2  2  2 
sin  if»  +  cos  ♦  sin  0 
to  to  to 

sin  5^/cos 

sin  a 

m 

-  sin  3a  cos  x. 

-  tan  *  /cos  0_ 

01  01 


Misslla  Angles  -  Earth  Angles 
0  -  9  -  sin”1  [sin  ylcon  i  ] 

10  ft  tfl 

i|i  «  sin”1  [sin  t|i  /cos  y  1  +  <5 

t0  ft  6 

Angular  Identities 
Earth  Angles 

tan  9  -  tan  x  cos  $ 

Missile  Velocities  -  Mlssila  Angles 

v  -  -  V  sin  * 

m 

w  ■  V  cos  sin  0 

m  m 

U  -  V  COS  *  cos  0 

Tm  m 
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Miaaile 


1 

ANGLE 

% 


Velocity  Ralatloaehlp 


COSINE 


TANGENT 


/  2  ,  2 
u  +  v 


07272 

U  +  V  +  w 


w 

u 

w 

u 

V 

u 
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Trajectory  Angle  -  Velocity  Relationshipa 


ANGLE 

SINE 

COSINE 

.2  .2  * 

Y. 

z 

(i  +  y  ) 

• 

V 

V 

5 

e 

y  _  . 

• 

X 

C 

2  2T 

2  IT 

(i  +  y  ) 

(i  +  y  ) 

TANGENT 


2  2H 
(x  +  y  ) 


x 
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APPENDIX  B 

6  DOF  COEFFICIENT  INPUT  LINE  NUMBERS 

AND 

ANALYSIS  INDEX  NUMBERS 
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APPENDIX  B 


6  DOF  Coefficient  Input  Line  Numbers 

and 

Analysis  Index  Numbers 


Each  coefficient  in  the  6  DOF  programs  has  an  associated  input  line  number  and 
analysis  index  number.  The  input  line  number  is  used  to  specify  initial 
conditions  at  fit  flags.  These  numbers  are  always  sequential  for  the  table 
currently  being  displayed.  Since  initial  fit  options  are  different  from 
subsequent  fit  options,  there  will  be  different  menus  for  each.  The  same  is 
true  for  fixed  plane  and  body  fixed  coordinate  systems. 

Within  the  analysis  program  there  is  one  vector  of  variables  for  all  the 
coefficients.  These  numbers  are  used  within  the  6  DOF  program  for  specifying 
parameters  and  fit  flags.  A  different  vector  exists  for  fixed  plane  and  body 
fixed  coordinate  systems. 

The  following  pages  will  define  the  parameters  available  within  the  6  DOF 
program  and  the  associated  line  and  index  numbers: 
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1.  FIXED  PLANE  LINE  NUMBERS 


* 


i 


* 


For  each  shot  in  the  analysis  (up  to  5),  a  menu  o£  aerodynamic 
coefficients  will  be  displayed.  For  the  last  shot  in  the  fit,  a  menu  of 
physical  parameters  is  displayed.  For  each  shot  fit,  menus  for  both 
aerodynamic  coefficients  and  physical  parameters  will  be  displayed. 

The  aerodynamic  coefficient  line  numbers  are: 


1  -  Cma 

2  -  Chkxg 

5  -  Cmqc*2 

6  -  Cnpa 

9  -  Cnsm 

10  -  Cma  mach 

13  -  Cxa^ 

14  -  CNa 

17  -  CYpa 

18  -  CYpa^ 

21  -  Cmq2-N 

22  -  CNa  mach 

25  -  Clp 

26  -  Clp  mach 

29  -  CYy«3 

30  -  Clyo^ 

33  —  Onttya 

34  -  Cnya 

3  -  Cmag 

4  -  Cmq 

7  -  Cnpotj 

8  -  Cnpa^ 

11  -  Cx 

12  -  Cx«2 

15  -  CNaj 

16  -  CNag 

19  -  Cx  mach 

20  -  Cmqa4 

23  -  Cnr2-N 

24  -  Clp 

27  -  CXyo2 

28  -  CZyoj 

31  -  Gnya^ 

32  -  Cnyotj 

The  physical  parameter  line  numbers 


1  -  Pitch  [THETA] 

2 

2  -  Yaw  [PSI] 

4 

5  -  Travel  [X] 

6 

7  -  Drift  [Y] 

8 

9  -  Height  [z] 

10 

11  -  Roll  [PHI] 

12 

13  -  Cma  unique 

14 

15  -  Cl  delta 

16 

17  -  CYyO 

18 

19  -  CnyO 

20 

21  -  Clp  unique 

are: 

-  Pitch  rate 

-  Yaw  rate 

-  Velocity 

-  Drift  rate 

-  Vertical  rate 

-  Roll  rate 

-  Inertia  ratio  [Ix/Iy] 

-  CZyO 

-  CmyO 

-  CXO  unique 
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2.  BODY  FIXED  LINE  NUMBERS 


For  each  shot  in  the  analysis  (19  to  5) ,  a  menu  of  aerodynamic 
coefficients  will  be  displayed.  For  the  last  shot  in  the  fit,  a  menu  of 
physical  parameters  is  displayed.  For  each  shot  fit,  menus  for  both 
aerodynamic  coefficients  and  physical  parameters  will  be  displayed. 

1  -  Craa  2  -  Cma-j  3  -  CN0  4  -  CN0^ 

5  -  Drag  6  -  Qraqc^  7  -  Cnr  8  -  Cnr2 


9  —  Crao^ 

10  -  CZotj 

11  -  Criya 

12  -  Cntya 

13  -  Cma-bar 

14  -  Cmctj-bar 

15  -  Cmq-bar 

16  -  Croqoij-bar 

17  -  Cnpo-bar 

18  -  CxO 

19  -  Cx«2 

20  -  Cx02 

21  -  Cx  roach 

22  -  CY0 

23  -  CY03 

24  -  CZa 

25  -  CZa 3 

26  -  CZya 

27  -  CYya 

28  -  CXa 

29  -  CYpa-bar 

30  -  CXo^-bar 

31  -  CNa-bar 

32  -  CNa^-bar 

33  -  Clp 

34  -  Clya 

35  -  C10 

36  -  Ixy 

37  -  Iy/Iz 

38  -  CmO 

39  -  CNO 

40  -  CYO 

41  -  CZO 

42  -  Cl-delta 

43  -  Cma-unique 

The  physical  parameter  line  numbers 

are: 

1  -  Pitch  [THETA] 

2  -  Pitch  rate 

2  -  Yaw  [PSI] 

4  -  Yaw  rate 

5  -  Travel  [X] 

6  -  Velocity 

7  -  Drift  [Y] 

8  -  Drift  rate 

9  -  Height  [z] 

10  -  Vertical  rate 

11  -  Roll  [PHI] 

12  -  Roll  rate 

13  -  Inertial  cross  product  [Ixy] 

14  -  Inertia  form  factor  [Ix/Iy] 

15  -  CmO 

16  -  CnO 

17  -  CYO 

18  -  CzO 

19  -  Cl  delta 

20  -  Cma-unique 

3.  FIXED  PLANE  INDEX  NUMBERS 
For  fixed  plane  equations  of 

are: 


1  -  Cma 

2  - 

5  -  Cmqo^ 

6  -  Cnpa 

9  -  Cnpctj 

10  -  Cma  mach 

17  -  Cx 

18  -  Cxotj 

21  -  CNo^ 

22  -  CNog 

25  -  Cx  mach 

26  -  Cmqa4 

29  -  Cnr2-N 

32  -  Clp 

35  -  CXyo^ 

36  -  Czya^ 

39  -  Cmyoc^ 

40  -  Cnyaj 

41  -  Pitch  [THETA] 

43  -  Yaw  [PSI] 

45  -  Travel  [X] 

47  -  Horizontal  motion  [Y] 

49  -  Vertical  motion  [Z] 

51  -  Roll  [PHI] 

53  -  Oma  unique 
55  -  Cl  delta 
57  -  CNdB 
59  -  CmdfiB' 

61  -  Clp  unique 
63  -  Yaw  rate  [RJ 
65  -  Hor.  rate  (missile)  [V] 
67  -  Roll  rate  [P] 


motion,  the  index  numbers  in  the  analysis 


3  -  Qmag 

4  -  Cmq 

7  -  Cnpa^ 

8  -  Cnpag 

11  -  Qmyog 

12  -  Cnyag 

19  -  Cxa4 

20  -  CNa 

23  -  CYpa 

24  -  CYpo^ 

27  -  Cmq2-N 

28  -  CNa  mach 

33  -  Clpo^ 

34  -  Clp  mach 

37  -  Cyyo^ 

38  -  Clyot2 

42  -  Pitch  rate 

44  -  Yaw  rate 

46  -  Velocity 

48  -  Horizontal  rate 

50  -  Vertical  rate 

52  -  Roll  rate 

54  -  Inertia  ratio  [Ix/Iy] 

56  -  CNdA 
58  -  CmdA 
60  -  CXO  unique 
62  -  Pitch  rate  [Q] 

64  -  Velocity  (missile)  [U] 

66  -  Vertical  rate  (missile)  [W] 
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4.  BODY  FIXED  INDEX  NUMBERS 


For  body  fixed  equations  of  motion,  the  index  numbers  in  the  analysis  are: 


1  -  Qua 

2  —  CmCTj 

3  -  CN3 

4 

-  CNPj 

5  -  Cmq 

6  -  Crnqc^ 

7  -  Cnr 

8 

*  Cnr2 

9  -  Qr  c  '2 

10  -  Czoj 

11  -  Cmyotj 

12 

-  Cnro^ 

13  -  Cma-bar 

14  -  Cmctj-bar 

15  -  Cmq-bar 

16 

-  Crtqoj-bat 

17  -  Cnpot-bar 

18  -  CxO 

19  -  0X0*2 

20 

- 

21  -  Cx  mach 

22  -  CY(3 

23  -  CY63 

24 

-  CZa 

25  -  CZotj 

26  -  CZyoc^ 

27  -  CYyo^ 

28 

-  CXa 

29  -  CYpa-bar 

30  -  CXot2“bar 

31  -  CNot-bar 

32 

-  CNotj-bar 

33  -  Clp 

34  -  Clr«2 

35  -  C10 

41  -  Pitch  [THETA] 

42  -  Pitch 

rate 

43  -  Yaw  [PSI] 

45  -  Travel  [X] 

47  -  Horizontal  motion  [Y] 

49  -  Vertical  motion  [Z] 

51  -  Roll  [PHI] 

53  -  Inertial  cross  product  [Ixy] 
55  -  CMO-trim 
57  -  CYO-trim 
59  -  Cl-delta 
61  -  Not  used 


44  -  Yaw  rate 

46  -  Velocity 

48  -  Horizontal  rate 

50  -  Vertical  rate 

52  -  Roll  rate 

54  -  Axial  Inertia  [lx] 

56  -  CNO-trim 

58  -  CZ 0-trim 

60  -  Cma-bar 

62  -  Pitch  rate  [Q] 


63  -  Yaw  rate  [R]  64  -  Velocity  (missile)  [U] 
65  -  Drift  rate  (missile)  [V]  66  -  Drop  rate  (missile)  [W] 
67  -  Roll  rate  [P] 


138 


